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Airborne  aerosols  contaminate  satellite  imagery,  making  it  difficult  for  analysts  to 
characterize  surface  targets.  The  amount  of  contamination  is  due  to  aerosol  loading,  best 
quantified  by  aerosol  optical  depth  (AOD).  Level  2  AERONET  data  provides  “ground 
truth”  AOD  measurements.  Unfortunately,  targets  are  not  likely  to  be  near  an 
AERONET  station.  This  thesis  examines  two  methods  of  quantifying  AOD  in  lieu  of 
AERONET  data:  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS)  AOD 
retrievals  and  Navy  Aerosol  Analysis  and  Prediction  System  (NAAPS)  AOD  forecasts 
(with  MODIS  data  assimilation).  Over  the  Persian  Gulf  and  Japan,  MODIS  accurately 
retrieved  AOD  at  0.55pm  wavelength.  MODIS  may  have  accurately  retrieved  AOD  over 
Korea,  but  uncertainty  in  the  retrieval  was  high  due  to  a  small  sample  size.  Over  West 
Africa,  MODIS  underestimated  AOD.  Over  the  Persian  Gulf  and  Japan,  NAAPS 
accurately  estimated  AOD  at  0.55pm.  Over  Korea,  the  model  displayed  mixed  results. 
Finally,  over  West  Africa,  NAAPS  tended  to  underestimate  AOD. 
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I.  INTRODUCTION 


A  major  problem  for  intelligence  agencies  that  analyze  satellite  images  is  the 
interaction  of  aerosol  particles  with  electro-optical  propagation.  If  atmospheric  aerosols 
obscure  a  target  scene,  image  analysts  will  likely  face  difficulty  in  characterizing  the 
target.  As  a  result,  depending  on  the  aerosol  contamination,  satellite  images  may  need 
some  correction  in  order  to  be  useful  to  the  intelligence  community. 

An  example  satellite  image  with  aerosol  contamination  is  shown  in  Figure  1.  The 
figure  displays  an  RGB  satellite  image  of  Iraq  taken  by  the  NOAA-17  satellite  on  26 
March  2003  (0718UTC).  A  large  dust  storm  covers  much  of  southwestern  Iraq,  making 
target  detection  and  target  characterization  difficult  in  this  region.  The  degree  of 
difficulty  in  target  detection  and  characterization  depends  on  the  aerosol  optical  depth. 


Figure  1.  NOAA-17  RGB  Satellite  Image  from  0718UTC,  26  March  2003.  A  large 
dust  storm  is  present  over  southwestern  Iraq  (From  NOAA  News  Online, 

Story  1120,2009) 


1 


The  analyst  must  be  able  to  correct  for  any  obscurations  covering  the  target,  in 
this  case  aerosols,  in  order  to  properly  detect  and  characterize  a  target.  Previous  studies 
by  have  used  Aerosol  Robotic  Network  (AERONET)  aerosol  optical  depth  (AOD)  data 
as  “ground-truth”,  indicating  that  AERONET  sites  reliably  provide  AOD  data  (Zhang  et 
ah,  2006,  2008).  Unfortunately,  targets  are  not  likely  to  be  located  near  an  AERONET 
site.  Previous  research  by  Misra  et  al.  (2008),  Chu  et  al.  (2002),  and  Zhang  et  al.  (2006) 
shows  that  polar  orbiting  satellite  instruments  such  as  the  Moderate  Resolution  Imaging 
Spectroradiometer  (MODIS)  can  also  provide  high  quality  AOD  data  for  a  target. 
However,  satellite -retrieved  AOD  data  is  only  available  at  the  times  that  the  satellite 
sensor  passes  over  the  target,  which  for  MODIS  is  only  twice  daily  during  daylight  hours. 
That  leaves  most  targets  with  22  hours  per  day  for  which  no  AOD  data  is  available.  The 
Navy  Aerosol  Analysis  and  Prediction  System  (NAAPS),  an  aerosol  model,  provides 
another  assessment  of  AOD  for  targets,  which  may  be  useful  when  no  other  AOD  data  is 
available. 

This  thesis  will  first  establish  an  AERONET  AOD  climatology  using  MODIS- 
retrieved  and  NAAPS-modeled  AOD  values  and  compare  these  techniques  against 
AERONET  AOD  measurements.  The  AERONET  AOD  climatology  will  show  the 
annual  and  seasonal  variability  for  aerosol  loading  for  each  of  the  four  basins:  Persian 
Gulf,  Korea,  Japan,  and  West  Africa.  The  comparison  of  these  techniques  against 
AERONET  AOD  measurements  will  test  if  each  technique  can  be  useful  as  a  proxy  for 
AERONET  AOD  measurements. 

Finally,  this  thesis  describes  a  radiance-based  target  detection  matrix.  To 
properly  use  this  matrix,  an  image  analyst  needs  information  about  the  target  and 
background  reflectance  (as  a  function  of  wavelength),  and  also  the  AOD.  The  radiance- 
based  target  detection  matrix  will  help  an  image  analyst  detennine  whether  or  not  a  target 
can  be  detected  and  characterized,  given  an  AOD  and  target/background  reflectance  as 
functions  of  wavelength. 
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II.  BACKGROUND 


A.  AEROSOLS,  SCATTERING,  AND  THE  VISIBLE  SPECTRUM 

The  focus  of  this  thesis  is  visible  solar  radiation  reflected  from  the  Earth’s 
surface.  Before  reaching  the  satellite,  reflected  solar  radiation  can  undergo  (in  many 
directions)  three  processes  within  the  atmosphere:  the  radiation  can  be  absorbed, 
scattered  by  gas  molecules  and  particles  (e.g.  aerosols,  cloud  droplets),  and/or  transmitted 
directly  to  the  satellite.  Kidder  and  Vonder  Haar  (1999)  noted  that  the  amount  of 
reflected  solar  radiation  by  a  particle  is  a  function  of:  particle  shape,  particle  size, 
particle  index  of  refraction,  wavelength  of  radiation,  and  viewing  geometry. 
Additionally,  in  1908,  Gustav  Mie  showed  that  the  scattered  radiation  is  a  function  only 
of  viewing  angle,  index  of  refraction,  and  the  size  parameter  (Kidder  1999).  The  size 
parameter  is  defined  as  (equation  3.49  from  Kidder  and  Vonder  Haar): 
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where  jis  the  size  parameter,  r  is  the  radius  of  the  scattering  particle,  and  2  is  the 
wavelength  of  the  radiation.  The  size  parameter  can  then  be  used  to  divide  scattering 
regimes  into  three  types:  Rayleigh,  Mie,  and  Geometric.  Mie  scattering  is  the  focus  of 
this  thesis  because  the  study  of  aerosols  (smoke,  dust,  haze)  using  the  visible  spectrum 
falls  within  the  Mie  regime. 

Aerosol  optical  depth  (AOD)  describes  the  reduction  of  photons  transmitted 
through  the  atmosphere  due  to  interactions  (i.e.,  scattering,  absorption)  with  aerosols. 
AOD  is  a  dimensionless  quantity,  and  is  a  function  of  wavelength.  Kidder  and  Vonder 
Haar  (1999)  define  optical  depth  as  (equation  3.28  from  Kidder  and  Vonder  Haar): 

zl 

Sx(z\,zi)  =  J  cre(/i  ,z)dz 

Zl 

where  8x  is  the  AOD  for  the  vertical  layer  extending  from  zi  to  zi .  The  right  side  of 

Equation  2  is  the  vertically  integrated  volume  extinction  coefficient  ( cr<? )  as  a  function  of 
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wavelength  for  the  layer.  An  increase  in  volume  extinction  coefficient  is  indicative  of 
increased  scattering  and/or  absorption  by  aerosol  particles.  Thus,  AOD  increases  as 
scattering  and/or  absorption  increases.  An  increase  in  scattering/absorption  occurs  when 
more  aerosol  particles  are  in  the  layer,  thereby  increasing  the  probability  of  photon 
interaction  with  aerosol  particles. 

Figure  2  illustrates  the  top  of  the  atmosphere  radiance  problem.  The  top  of  the 
atmosphere  solar  irradiance  is  denoted  as  Eo,  and  the  solar  irradiance  at  the  Earth’s 
surface  is  represented  by  Es.  The  top  of  the  atmosphere  irradiance  is  greater  than  the 
amount  of  solar  irradiance  at  the  Earth’s  surface  due  to  scattering  and/or  absorption  of  the 
irradiance  as  the  photons  pass  through  the  atmosphere.  The  radiance  reflected  by  the 
target  is  Lt,  and  the  radiance  reflected  by  the  background  is  Lb.  The  amount  of  radiance 
reflected  by  the  target  and  background  is  a  function  of  wavelength  and  surface 
characteristics  (i.e.,  surface  color). 

Photons  reflected  by  the  surface  can  be  absorbed  or  scattered  in  many  directions 
as  they  pass  through  the  atmosphere  on  the  way  to  the  satellite.  In  order  for  a  satellite- 
borne  sensor  to  distinguish  a  target  from  its  background,  the  difference  between  the  target 
and  background  top  of  the  atmosphere  radiance  (A Lt)  must  be  greater  than  the  sensor 
radiance  sensitivity  ( A L? ).  The  difference  between  the  target  and  background  top  of  the 
atmosphere  radiance  is  defined  by  the  following  equation  (Durkee  2009): 

alt=  Lrr  LbT  =  (Lt'  Lb)e  5 

The  top  of  the  atmosphere  radiance  from  the  target,  LtT,  is  the  transmitted  radiance 
reflected  by  the  target  surface  and  the  path  radiance.  The  top  of  the  atmosphere  radiance 
from  the  background,  LbT,  is  the  transmitted  radiance  reflected  by  the  background  surface 
and  the  path  radiance.  Aerosol  concentration  (i.e.,  AOD)  is  important  because  high 
aerosol  concentrations  over  a  target  can  increase  the  amount  of  scattering  and/or 
absorption  of  photons  reflected  by  the  target,  thereby  reducing  the  contrast  between  the 
target  and  its  background.  Reduced  target/background  contrast  lowers  the  probability 
that  a  satellite  sensor  will  distinguish  between  the  target  and  background. 
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Figure  2.  Radiance  received  by  a  satellite  at  the  top  of  the  atmosphere,  with  top  of  the 
atmosphere  irradiance  (E0),  solar  irradiance  at  the  Earth’s  surface  (Es),  optical 
depth  (8),  and  reflected  radiance  by  the  target  (Lt)  and  background  (Lb)  noted 
above.  (From  Durkee  February  2009) 

B.  TROPOSPHERIC  DUST  AND  SMOKE  CLIMATOLOGY 

Climatology  narratives  are  included  to  describe  seasonal  variations  in  weather 
patterns  for  the  Persian  Gulf,  East  Asia,  and  West  Africa.  The  following  narratives  are 
paraphrased  from  the  climatology  narratives  produced  by  the  Fourteenth  Weather 
Squadron  in  Asheville,  NC. 

1.  Persian  Gulf 

The  climate  of  the  Persian  Gulf,  specifically  that  of  Saudi  Arabia  and  the  United 
Arab  Emirates  (UAE),  has  characteristics  of  a  monsoon  climate,  because  the  area  falls 
under  the  influence  of  the  Near  Equatorial  Tradewind  Convergence  (NETWC)  zone.  The 
NETWC  boundary  stays  south  of  the  UAE,  so  the  prevailing  year-round  wind  flow  for 
this  region  is  north-northwesterly.  The  strength  of  the  flow  varies  depending  on  the  time 
of  year,  which  effectively  leads  to  two  seasons:  Winter  and  Summer. 
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The  Winter  season  covers  November  through  April.  During  Winter,  restrictions 
to  visibility  occur  due  to  the  presence  of  dust,  sand,  haze,  and  fog  (along  the  coast). 
Overall,  blowing  sand  and  dust  reduces  surface  visibilities  below  seven  miles 
approximately  three  times  per  month.  However,  if  the  flow  is  stronger  and  the  dust/sand 
is  lofted  higher  in  the  atmosphere,  visibilities  can  be  reduced  for  days  at  a  time,  which 
can  hinder  target  detection  for  satellite-borne  sensors. 

The  Summer  season  extends  from  May  through  October.  Flow  remains  north- 
northwesterly,  with  prevailing  surface  visibilities  between  four  and  seven  miles  in  dust, 
sand,  or  haze.  Stronger  flow  over  northern  Saudi  Arabia  can  loft  so  much  dust  into  the 
atmosphere  that  the  dust  remains  suspended  for  weeks  at  a  time,  further  hindering  target 
detection  for  satellite  sensors.  Often  accompanying  the  strong  northwesterly  flow  is  a 
low-level  Persian  Gulf  jet,  which  further  enhances  the  amount  of  dust  lofted  over  the 
western  Persian  Gulf. 

Shamal  winds  can  occur  anytime  during  the  year  but  are  especially  common 
during  the  Summer  months.  Shamal  winds  are  strong  northwesterly  wind  events  that 
originate  over  the  Tigris-Euphrates  lowlands  in  Iraq.  Shamal  events  often  carry  large 
amounts  of  dust  and  sand,  lasting  for  periods  of  one  to  five  days.  However,  some  long 
lived  Shamal  wind  storms  have  lasted  40  days  during  the  summer,  which  can  play  havoc 
with  attempts  to  image  targets  via  satellite  over  this  region.  Consequently,  Summer  is  the 
season  with  the  heaviest  aerosol  loading  over  the  Persian  Gulf. 

2.  East  Asia 

Korea  and  Japan  experience  four  distinct  seasons.  The  first  season  is  Winter, 
which  stretches  from  November  through  March.  During  Winter,  the  Asiatic  high  over 
China  builds  in  November,  reaches  peak  strength  in  January,  and  weakens  in  March,  with 
predominately  northwesterly  flow.  The  Asiatic  high  is  the  dominating  feature  over  East 
Asia  during  the  winter.  Consequently,  cold,  dry,  continental  polar  air  flows  southward 
from  China  over  Korea  and  Japan.  Clear  skies  and  little  restriction  to  visibility  are  the 
main  results  of  the  influence  of  the  Asiatic  high.  Occasionally,  however,  dust  stonns 
(Yellow  Wind)  from  China  and  Mongolia  reduce  visibility  over  Korea  and  Japan. 
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Yellow  Wind  events  occur  one  to  two  times  a  month  from  December  through  February 
and  two  to  three  times  per  month  during  March,  though  less  frequently  over  Japan.  Smog 
around  the  cities  also  reduces  surface  visibilities  during  the  morning.  Both  Yellow  Wind 
events  and  smog  hinder  surface  visibility  and  can  negatively  impact  the  ability  of  satellite 
sensors  to  detect  a  target  on  the  ground. 

Spring  only  lasts  from  April  through  May,  and  is  the  transition  between  the  clear, 
cold,  dry  Winter  season  to  the  cloudy,  wann,  wet  summer  season.  During  Spring, 
northwesterly  flow  begins  to  weaken  as  the  Asiatic  high  weakens  and  is  displaced 
northwards.  Also,  more  frontal  systems  move  over  East  Asia  and  the  flow  becomes 
westerly.  Consequently,  as  frontal  passages  become  more  common,  so  does  the 
frequency  of  Yellow  Wind  events,  especially  in  April,  which  hinders  surface  visibility 
and  target  detection.  Again,  as  in  winter,  smog  around  the  cities  also  reduces  surface 
visibilities  in  the  morning  hours. 

Summer  over  East  Asia  extends  from  June  until  mid-September.  Summer  is  the 
cloudy,  warm,  and  wet  season,  with  the  weather  largely  determined  by  the  position  of  the 
monsoon  frontal  boundary.  The  monsoon  frontal  boundary  separates  the  moist,  wann, 
maritime  tropical  air  mass  from  the  cooler,  dry,  continental  air  mass,  and  brings  with  it 
lots  of  precipitation,  cloud  cover,  and  fog.  Restrictions  to  visibility  during  the  Summer 
months  are  predominately  due  to  clouds,  precipitation,  and  fog  rather  than  smog  and  dust. 

Fall  covers  mid-September  to  November  and  is  the  transition  season  from  the 
cloudy,  warm,  wet  summer  season  to  the  clear,  cold,  dry  winter  season.  As  the  Asiatic 
high  re-establishes  itself,  the  flow  becomes  more  westerly  and  the  number  of  days  with 
clouds  and  precipitation  begins  to  decrease.  Fog  and  smog  are  the  main  restrictors  to 
surface  visibility  during  the  Fall  season. 

3.  West  Africa 

The  climate  of  Benin  and  Nigeria  is  governed  by  the  position  of  the  monsoon 
trough.  The  monsoon  trough  in  this  region  is  better  known  as  the  NETWC  zone.  The 
NETWC  separates  the  hot,  dry,  continental  air  mass  that  dominates  the  Sahara  from  the 
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warm,  moist,  maritime  air  mass  that  dominates  the  tropics.  As  the  NETWC  moves  north 
or  south,  so  do  these  air  masses,  giving  rise  to  two  seasons  for  Benin  and  Nigeria.  These 
seasons  are  the  dry  season  and  the  wet  season. 

The  dry  season  spans  November  through  March,  roughly  Northern  Hemisphere 
winter.  During  the  dry  season,  the  NETWC  moves  south  towards  the  coastline  of  these 
two  countries,  and  this  ushers  in  dry  Harmattan  winds.  Harmattan  winds  are  northerly 
winds  that  carry  Saharan  dust  and  sand  southward,  which  lowers  visibility  and  obscure 
targets  from  satellite  observation.  In  fact,  Harmattan  winds  can  reduce  visibilities  to  one- 
half  mile  up  to  altitudes  of  12,000  feet,  though  more  often  visibility  is  restricted  to  less 
than  three  miles  up  to  altitudes  of  8,000  feet.  Consequently,  the  dry  season  is  the  time  of 
year  that  aerosol  loading  is  heaviest  over  West  Africa. 

The  wet  season  spans  April  through  October,  roughly  Northern  Hemisphere 
summer.  During  the  wet  season,  the  NETWC  moves  north  of  Benin  and  Nigeria,  which 
allows  the  prevailing  flow  to  turn  south-southwesterly.  As  the  flow  becomes  more 
southerly,  wann,  moist,  maritime  air  is  brought  into  the  area.  Thus,  any  restrictions  to 
visibility  and  target  detections  are  primarily  due  to  clouds,  precipitation,  and/or  fog. 

In  addition  to  sand  and  dust,  biomass  burning  is  another  source  of  aerosols  in 
West  Africa  that  can  reduce  visibility  and  target  detection.  Biomass  burning  is  especially 
common  in  the  savanna  region,  where  land  is  burned  after  harvest,  usually  at  the  end  of 
the  rainy  season  (NASA  2009). 
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III.  DATA  AND  METHODOLOGY 


A.  AEROSOL  ROBOTIC  NETWORK  (AERONET) 

1.  Network  Description 

AERONET  is  a  network  of  ground-based  sun  photometers  established  by  NASA 
and  the  French  National  Center  for  Scientific  Research  (CRNS)  in  collaboration  with 
many  universities  and  agencies  worldwide.  The  purpose  of  this  network  is  to  provide 
standard,  continuous,  and  long-term  worldwide  measurements  of  aerosol  optical  depths, 
inversion  products,  and  aerosol  regimes.  The  AERONET  sun  photometers  measure  the 
intensity  of  direct  sunlight  and  directional  sky  brightness  throughout  nine  bands  from 
0.34pm  (UV)  to  1.64pm  (near  IR).  Data  are  then  relayed  by  satellite  to  the  NASA 
Goddard  Space  Flight  Center,  where  it  is  quality  controlled  and  processed  into  the 
products  listed  above.  For  this  thesis,  Level  2  data  are  used  as  ground  truth  for  eight 
AERONET  sites  throughout  the  West  Africa,  Persian  Gulf,  and  East  Asia  basins.  Level  2 
data  are  the  highest  quality  data,  as  a  result  of  being  both  cloud  screened  and  quality 
assured.  AERONET  AOD  data  has  an  uncertainty  of  0.01  to  0.02  in  AOD  (GFSC  2009). 

2.  Site  Selection 

The  first  consideration  in  choosing  AERONET  sites  for  this  study  is  the 
availability  of  Level  2  AERONET  data.  Data  prior  2006  are  not  available  for  the 
AERONET  sites  used  in  this  thesis.  The  2006-2007  period  provides  sufficient  data  for 
comparisons  with  NAAPS  and  MODIS.  For  the  Persian  Gulf  basin,  Level  2  AERONET 
AOD  data  are  available  on  80%  of  the  days  during  the  two-year  period  (2006-2007). 
Level  2  AERONET  AOD  data  are  available  only  on  30%  and  43%  of  the  days  for  the 
Korea  and  Japan  basins,  respectively.  For  the  West  Africa  basin,  data  are  available  on 
68%  of  the  days  during  the  two-year  period  (2006-2007). 
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The  second  consideration  in  choosing  AERONET  sites  for  the  thesis  is  the 
proximity  of  the  AERONET  sites  to  the  coast.  The  MODIS  AOD  data  used  in  this  study 
is  over-ocean  data.  In  order  to  have  valid  comparisons  between  the  AERONET  and 
MODIS  AOD  data,  AERONET  sites  need  to  be  close  to  the  coast.  Of  the  eight 
AERONET  sites  this  thesis  includes,  five  of  them  are  located  along  the  coast.  The 
coastal  AERONET  sites  are:  Dhabi,  UAE,  Dhadnah,  UAE,  Anmyon,  South  Korea, 
Osaka,  Japan,  and  Shirahama,  Japan.  The  inland  AERONET  sites  are:  Solar  Village, 
Saudi  Arabia,  Djougou,  Benin,  and  Ilorin,  Nigeria. 


The  locations  of  the  AERONET  sites  throughout  the  four  basins  are  identified  in 
Figures  3,  4  and  5.  Figure  3  displays  the  locations  of  AERONET  sites  in  the  Persian  Gulf 
basin.  The  Dhabi  and  Dhadnah  AERONET  sites  are  along  the  coast,  while  Solar  Village 
is  approximately  417  km  inland.  Figure  4  displays  the  locations  of  AERONET  sites  in 
East  Asia.  All  of  the  East  Asian  AERONET  sites  used  in  this  study  are  located  along  the 
coast.  Figure  5  displays  the  locations  of  the  AERONET  sites  in  the  West  Africa  basin. 
Both  of  the  West  African  AERONET  sites  are  considerably  inland  from  the  Gulf  of 
Guinea,  with  Ilorin  and  Djougou  approximately  217  km  and  400  km  inland,  respectively. 


A 

<  >) 


AMyftrV  - 
Irydaf- 

-J 


Ad  Dlr'iyah 

°  o 

Ar  Rlyad  o 


Al  Qatif. 
Sihat 


Ad  Dammam. 
Al  Khobar 


Hofuf. 

Mubarraz 


^  Map  |  Sateinte  |  Ten<iin  | 

Dhadnah,.  UAE 

Dhabi, _UAE- 

Solar  Village,  Saudi 
Arabia 


Figure  3.  Three  Persian  Gulf  AERONET  sites  used  in  this  study  (From  Google  Maps 

February  2009). 
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Figure  4.  Three  East  Asia  AERONET  sites  used  in  this  study  (From  Google  Maps 

February  2009). 
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3. 


AERONET  Wavelength  Variability 


The  NAAPS  and  MODIS  AOD  values  are  retrieved  at  0.55pm.  However,  no 
AOD  data  at  0.55)iun  are  available  at  any  of  AERONET  sites.  AERONET  AOD  data  are 
typically  available  at  0.34|um,  0.38pm,  0.44pm,  0.50pm,  0.675pm,  0.87pm,  1.02pm,  and 
1.675pm.  This  thesis  compares  the  NAAPS  and  MODIS  AOD  data  to  the  0.50pm, 
AERONET  AOD  data,  the  closest  available  ground-truth  wavelength  (note:  0.44pm 
wavelength  is  used  for  Djougou,  since  no  0.50pm  AOD  data  are  available).  Figure  61 
(Appendix  A)  displays  the  differences  in  AOD  distributions  by  wavelength  for  Solar 
Village.  Differences  in  AOD  values  between  the  wavelengths  are  small.  For  example,  at 
an  AOD  of  0.40,  the  cumulative  frequencies  at  0.50pm  and  0.675pm  are  72%  and  78%, 
respectively,  only  a  6%  difference.  Therefore,  this  study  will  use  AERONET  AOD 
measurements  at  0.50pm  as  a  proxy  for  AOD  measurements  at  0.55pm. 

B.  MODERATE  RESOLUTION  IMAGING  SPECTRORADIOMETER 

(MODIS) 

1.  MODIS  Data  Description 

MODIS  is  a  key  part  of  the  mission  of  the  NASA  EOS  program.  The  MODIS 
sensor  flies  on  two  satellites  (Aqua  and  Terra)  that  together  provide  global  coverage  over 
a  one  to  two  day  time  period.  Each  satellite  is  located  at  an  altitude  of  705km,  has  a  view 
swath  width  of  2,330km,  and  scans  the  Earth  across  36  spectral  bands,  ranging  from 
0.405  to  14.385pm.  Terra  follows  a  descending  orbital  path,  passing  over  the  Equator  at 
1030LT  (Local  Time),  while  Aqua  follows  an  ascending  orbital  path,  passing  over  the 
Equator  at  1330LT.  With  one  morning  and  one  afternoon  pass  over  the  same  area,  it  is 
possible  to  detect  changes  in  the  atmosphere  that  occur  during  the  day,  in  addition  to 
minimizing  the  effects  of  sun-glint  and  shadows  at  the  surface  (GFSC  2009). 

The  Naval  Research  Laboratory,  Monterey,  CA,  retrieves  aerosol  optical  depth 
using  the  MOD04/MYD04  algorithm  for  the  0.55pm  (green)  spectral  band  (Zhang  et  al. 
2006).  The  MOD04/MYD04  algorithm  retrieves  over-ocean  AOD  values  and  bins  the 
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AOD  values  into  1°  x  1°  latitude/longitude  grids,  and  does  not  retrieve  AOD  over  land 
(Zhang  et  al.  2008).  The  MODIS  AOD  grid  point  value  is  the  mean  of  the  values  in  the 
1°  x  1°  box.  This  thesis  uses  daily  mean  over-ocean  AOD  values  retrieved  by  NRLMRY 
for  2006-07.  Clouds  and  sun  glint  prevent  an  AOD  retrieval.  The  NRL  dataset, 
therefore,  does  not  provide  a  daily  mean  value  for  every  grid  point.  The  uncertainty 
associated  with  the  MODIS  derived  aerosol  optical  depths  is  0.03  ±  0.055  (Remer  et  al. 
2005). 


2.  MODIS  Climatology 

Histograms  of  MODIS  AOD  are  prepared  for  each  basin  for  the  two-year  period 
and  for  each  season.  Clouds  and  sun  glint  prevent  a  daily  mean  AOD  value  to  be 
available  at  the  grid  point  closest  to  an  AERONET  station.  Therefore,  AOD  data  from 
several  grid  points  near  each  AERONET  site  are  used  to  build  the  MODIS  AOD 
climatology  histograms  for  each  basin.  The  MODIS  histograms  will  be  compared  with 
the  AERONET  histograms  by  basin. 

3.  MODIS/ AERONET  Comparison  Methodology 

In  addition  to  comparing  MODIS  and  AERONET  basin  histograms,  point-to- 
point  MODIS/AERONET  comparisons  are  made.  The  first  step  is  to  make  sure  that  the 
MODIS  and  AERONET  data  are  matched  in  space.  The  MODIS  grid  point  that  is  closest 
to  the  AERONET  site  is  selected.  The  second  step  is  to  make  sure  that  the  MODIS  and 
AERONET  data  are  matched  in  time  for  each  site.  To  match  the  time  of  MODIS 
overpasses  (approximately  1030  and  1330L),  AERONET  measurements  from  1000L  to 
1400L  are  used  to  calculate  a  mean  AERONET  AOD  value.  The  daily  mean  MODIS 
AOD  value  is  then  compared  with  the  mean  AERONET  AOD  value  for  each  day  over 
the  two-year  period  to  produce  seasonal  and  annual  scatter  plots. 

There  are  days  when  AERONET  measurements  are  not  taken,  because  of  cloud 
cover,  malfunctioning  equipment,  etc.  Additionally,  there  are  some  days  for  which 
MODIS  data  is  unavailable.  If  either  AERONET  or  MODIS  data  is  unavailable  for  a 
specific  day,  that  day  is  not  included  in  this  study. 
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4. 


MODIS  Grid  Locations 


Figure  6  illustrates  the  MODIS  AOD  grid  points  in  the  Persian  Gulf  basin  which 
are  close  to  AERONET  sites  Dhabi,  UAE,  and  Dhadnah,  UAE.  MODIS  grid  points  are 
the  numbered.  The  13  grid  points  are  used  to  create  the  Persian  Gulf  basin  MODIS  AOD 
climatology  histograms.  For  the  point-to-point  MODIS/AERONET  comparisons, 
MODIS  AOD  data  from  Grid  Point  1  (red)  in  the  Persian  Gulf  is  compared  to  AERONET 
AOD  data  from  Dhabi  and  MODIS  AOD  data  from  Grid  Point  1  (green)  in  the  Gulf  of 
Oman  are  compared  to  AERONET  AOD  data  from  Dhadnah.  Data  from  the  Solar 
Village,  Saudi  Arabia  AERONET  site  are  not  included  because  the  AERONET  site  is  too 
far  inland  (417  km)  to  have  a  valid  comparison  with  over-ocean  MODIS  data. 

Figure  7  illustrates  the  MODIS  grid  points  for  the  Korea  and  Japan  basins  which 
are  close  to  Anmyon,  S.  Korea  and  Shirahama,  Japan  AERONET  sites.  The  nine 
numbered  grid  points  (red)  are  used  to  create  the  Korea  basin  MODIS  AOD  climatology 
histograms,  while  the  six  numbered  grid  points  (green)  are  used  to  create  the  Japan  basin 
histograms.  For  the  point-to-point  MODIS/AERONET  comparisons,  MODIS  data  from 
Grid  Point  1  (red)  are  compared  to  AERONET  data  from  Anmyon  to  produce  the  Korea 
scatter  plots.  Similarly,  MODIS  data  from  Grid  Point  1  (green)  are  compared  to 
AERONET  data  from  Shirahama  to  produce  the  Japan  scatter  plots.  Data  from  the  Osaka 
(Japan)  AERONET  site  are  not  included  due  to  its  small  dataset  for  2006-2007. 

Figure  8  illustrates  the  MODIS  grid  points  for  the  West  Africa  basin  which  are 
closest  to  the  Ilorin,  Nigeria  AERONET  site.  The  10  numbered  pixels  are  used  to  create 
the  West  Africa  basin  MODIS  AOD  climatology  histograms.  MODIS  AOD  data  from 
grid  point  4  are  compared  to  AERONET  AOD  data  from  Ilorin  to  produce  the  West 
Africa  scatter  plots.  Data  from  the  Djougou,  Benin  AERONET  site  are  not  included 
because  the  AERONET  site  is  far  too  inland  (400  km)  to  make  any  valid  comparisons 
with  over-ocean  MODIS  data. 


14 


<0> 


6 

Ad  Oawftah 


■*><5 


Dhadnah,  UAE 
Dhabi,  UAE 


1  0 


United  Arab 
Emirates 


o  Masqat 

o  0  O 

AlAub* 


MtrlaUa 

■M'/'.j'V 


Figure  6.  MODIS  grid  points  and  AERONET  sites  used  for  the  Persian  Gulf  basin, 
MODIS/AERONET  comparison.  The  numbers  denote  MODIS  grid  points, 
and  the  AERONET  sites  are  denoted  by  filled  circles.  The  six  grid  points  in 
red  are  used  to  represent  the  over-water  region  near  Dhabi  (red  circle).  The 
seven  grid  points  in  green  are  used  to  represent  the  over-water  region  near 
Dhadnah  (green  circle)  (From  Google  Maps  February  2009). 


Figure  7.  Same  as  Figure  6,  but  for  the  Korea  (red)  and  Japan  (green)  basins.  The 
AERONET  sites  are  denoted  by  filled  circles  (From  Google  Maps  February 

2009). 
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Figure  8.  Same  as  Figure  6,  but  for  West  Africa.  The  AERONET  site  (Ilorin)  is 
denoted  by  the  red  circle  (From  Google  Maps  February  2009) 

C.  NAVY  AEROSOL  ANALYSIS  AND  PREDICTION  SYSTEM  (NAAPS) 

1.  NAAPS  AOD  Description 

NAAPS  is  the  model  developed  by  the  Naval  Research  Laboratory  (NRL)  in 
Monterey,  CA,  and  is  used  to  predict  the  distribution  of  tropospheric  aerosols.  NAAPS 
uses  the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  to  provide 
global  meteorological  fields.  Twice  daily,  NAAPS  produces  analyses  and  forecasts  of 
aerosol  concentrations  at  six  hour  intervals  for  a  five-day  forecast  period  on  1°  x  1° 
latitude/longitude  grids  (NAAPS  2009  and  Zhang  et  al.  2008).  Aerosol  optical  depth  is 
not  directly  calculated  by  NAAPS.  Aerosol  optical  depth  is  produced  by  the  Forecast  of 
Aerosol  Radiative  Optical  Properties  (FAROP)  optical  model.  NAAPS  uses  the  FAROP 
optical  model  as  a  post-processor,  which  takes  the  aerosol  concentration  data  from 
NAAPS  model  runs  and  computes  optical  properties  of  aerosols,  i.e.,  aerosol  optical 
depth. 
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This  thesis  uses  the  AOD  product  produced  from  NAAPS  with  assimilated  over¬ 
ocean  MODIS  AOD  data  using  the  Navy  Atmospheric  Variational  Data  Assimilation 
System  (NAVDAS).  NAVDAS  assimilates  over-ocean  MODIS  AOD  data  at  0.55pm  by 
first  mapping  the  MODIS  data  to  1°  x  1°  grids.  To  prevent  AOD  errors  introduced  by  the 
presence  of  clouds,  two  quality  control  procedures  are  applied  to  the  MODIS  AOD  data. 
First,  when  the  cloud  fraction  is  larger  than  30%  for  a  grid  box,  the  grid  point  is  excluded 
from  the  dataset.  Second,  bins  with  less  than  five  data  entries  and  bins  isolated  from 
continuous  aerosol  features  are  removed  (Zhang  et  al.  2008). 

The  NAAPS  AOD  product  used  in  this  thesis  is  an  analysis  on  a  1°  x  1°  grid  at 
00,  06,  12  and  18  UTC  for  each  day  in  2006-07.  NAAPS  AOD  values  are  available  for 
over-ocean  points  only.  There  is  a  NAAPS  AOD  grid  point  within  30km  of  AERONET 
sites  Dhabi,  UAE,  Anmyon,  South  Korea,  Osaka  and  Shirahama,  Japan.  The  closest  grid 
point  to  AERONET  site  Dhadnah,  UAE,  is  1.0°  longitude  due  east  of  the  site.  The  Ilorin, 
Nigeria  AERONET  site  is  inland  from  the  coast.  The  over-ocean  NAAPS  AOD  grid 
point  used  to  compare  with  Ilorin  is  2.8°  latitude  due  south  of  the  AERONET  site.  The 
NAAPS  AOD  product  provides  an  AOD  value  for  each  over-ocean  grid  point  at  each 
analysis  time.  The  NAAPS  AOD  grid  point  closest  to  the  six  AERONET  sites  will  be 
used  in  annual  and  seasonal  histograms  and  also  in  point-to-point  comparisons  with 
AERONET.  Two  AERONET  sites,  Djougou,  Benin,  and  Solar  Village,  Saudi  Arabia, 
are  not  included  in  the  NAAPS/AERONET  comparisons  because  both  sites  are  too  far 
inland. 


2.  NAAPS/AERONET  Comparison  Methodology 

Histograms  of  NAAPS  AOD  are  prepared  for  each  basin  for  the  two-year  period 
and  for  each  season.  In  addition  to  comparing  NAAPS  and  AERONET  basin  histograms, 
point-to-point  NAAPS/AERONET  comparisons  are  made.  For  each  AERONET  site,  a 
daily  mean  AERONET  AOD  value  is  computed.  Using  the  closest  NAAPS  grid  point  to 
each  AERONET  site,  the  daily  mean  NAAPS  AOD  value  is  also  computed.  On  days 
when  both  AERONET  and  NAAPS  daily  mean  values  are  available,  a  point-to-point 
comparison  is  perfonned. 
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D.  RADIANCE-BASED  TARGET  DETECTION 


1.  MODTRAN  ™  Model  Runs 

MODTRAN™  is  the  U.S.  Air  Force  moderate  spectral  resolution  radiative 
transport  model  for  wavelengths  in  the  ultraviolet  to  thermal  IR.  In  this  thesis, 
MODTRAN™  is  used  to  compute  the  top-of-the-atmosphere  (TOA)  radiance  for  an  input 
vertical  profile  of  atmospheric  temperature,  moisture  and  dust  aerosols  typical  of  a  desert 
environment.  The  surface  reflectance  is  an  input  to  MODTRAN™.  To  simulate  ground- 
based  targets  of  different  reflectance  characteristics,  multiple  MODTRAN™  runs  are 
made,  each  with  a  different  surface  reflectance  value  and  the  same  amount  of  dust 
aerosols.  To  simulate  the  radiative  effects  of  different  aerosol  loading,  multiple 
MODTRAN™  runs  are  made,  each  with  a  different  amount  of  dust  aerosols  and  the  same 
surface  reflectance.  In  all,  the  MODTRAN™  model  is  run  260  times  with  different 
surface  reflectance  values  and  dust  aerosol  loading.  The  model  runs  simulate  the  TOA 
radiance  expected  to  be  detected  by  a  satellite  sensor  for  AOD  values  ranging  from  0.0- 
3.5  and  surface  reflectance  values  ranging  from  0-99%. 

2.  MODTRAN™  Radiance  Comparison  with  MODIS  Imagery 

This  study  compares  MODTRAN™  TOA  radiance  estimates  with  TOA  radiance 
in  MODIS  imagery.  The  MODIS  radiance  values  are  treated  as  ground  truth  for  the 
comparison.  The  purpose  of  the  MODTRAN/MODIS  comparison  is  to  establish  the 
reliability  MODTRAN™  TOA  radiance  values.  Then,  MODTRAN™  can  be  used  to 
create  a  radiance-based  target  detection  matrix.  Results  of  this  analysis  are  discussed  in 
Chapter  IV,  Section  D. 


18 


IV.  RESULTS 


A.  AERONET  CLIMATOLOGY 

1.  Basin  Composite 

This  thesis  establishes  a  basin  composite  AERONET  climatology.  The  purpose  is 
to  determine  how  the  different  geographic  regions  rank  in  terms  of  aerosol  optical  depth. 
First,  this  study  needs  to  detennine  if  there  are  any  significant  differences  in  the  basin 
composite  AERONET  AOD  distributions  for  2006  and  2007.  Figure  9  is  the  one-year 
AERONET  AOD  (0.50pm)  distribution  for  2006  using  data  from  the  Persian  Gulf, 
Korea,  Japan,  and  West  Africa  basins.  The  figure  contains  both  a  histogram  of  AOD 
values  and  cumulative  frequency  distribution  corresponding  to  the  histogram.  Half  of  the 
AOD  distribution  lies  between  values  of  0.21  and  0.56,  with  the  median  of  the 
distribution  at  0.35. 

Figure  10  is  the  one -year  AERONET  0.50pm  wavelength  AOD  distribution  for 
2007  using  data  from  the  four  basins  of  interest.  The  figure  contains  both  a  histogram  of 
AOD  values  and  cumulative  frequency  distribution  corresponding  to  the  histogram.  Half 
of  the  AOD  distribution  lies  between  values  of  0.21  and  0.59,  with  the  median  of  the 
distribution  at  0.39. 
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One-year  distribution  of  basin  composite  AERONET  AOD  (0.50pm)  values 
for  the  Persian  Gulf,  Korea,  Japan,  and  West  Africa  basins,  for  2006.  The 
histogram  of  AOD  values  (blue  bars)  and  cumulative  frequency  distribution 
(red  curve)  are  denoted.  The  sample  size  is  listed  in  the  upper  right  corner. 

The  vertical  lines  (black),  from  left  to  right,  denote  the  quartiles  of  the 
distribution,  25%  (dashed),  50%  (solid),  75%  (dashed). 
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Figure  10.  Same  as  Figure  9,  but  for  2007. 
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Figure  1 1  illustrates  the  differences  between  the  AERONET  AOD  distributions 
for  2006  and  2007.  The  cumulative  frequency  curve  for  each  individual  year  is  shown. 
The  cumulative  frequency  curve  for  2006  closely  mirrors  the  cumulative  frequency  curve 
from  2007,  indicating  no  significant  differences  in  the  distributions  for  each  year. 
Therefore,  the  2006  and  2007  AERONET  AOD  datasets  can  be  safely  combined  to 
produce  a  basin  composite  AOD  climatology. 


Figure  1 1 .  Cumulative  frequency  curves  for  the  basin  composite  climatology  for  2006- 
2007.  The  cumulative  frequency  curves  for  2006  (blue)  and  2007  (green)  are 
are  computed  from  the  distribution  for  each  year.  The  combined  cumulative 
frequency  curve  (red)  is  computed  from  the  combined  dataset  from  both  years. 


Figure  12  is  the  two-year  composite  AERONET  0.50pm  AOD  distribution  for 
2006-2007  using  AERONET  AOD  data  from  the  four  basins.  The  figure  contains  both  a 
histogram  of  AOD  values  and  cumulative  frequency  distribution  corresponding  to  the 
histogram.  Half  of  the  AOD  distribution  lies  between  AOD  values  of  0.21  and  0.59,  with 
the  median  of  the  distribution  at  0.38.  In  6%  of  the  cases,  AOD  values  are  greater  than 
1.0. 
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To  demonstrate  why  the  2006-2007  distribution  could  be  useful  for  an  image 
analyst,  assume  an  AOD  equal  to  or  greater  than  0.50  makes  an  image  uncorrectable. 
Figure  12,  therefore,  suggests  that  aerosol  loading  may  be  light  enough  to  allow  for  scene 
correction  in  approximately  65%  of  the  cases. 
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Figure  12.  Two-year  distribution  of  basin  composite  AERONET  AOD  (0.50pm)  values 
for  the  Persian  Gulf,  Korea,  Japan,  and  West  Africa  basins,  for  2006-2007. 

The  histogram  of  AOD  values  (blue  bars)  and  cumulative  frequency 
distribution  (red  curve)  are  denoted.  The  sample  size  is  listed  in  the  upper 
right  corner.  The  vertical  lines  (black),  from  left  to  right,  denote  the  quartiles 
of  the  distribution,  25%  (dashed),  50%  (solid),  75%  (dashed). 

In  terms  of  relative  contributions  to  the  composite  dataset,  the  Persian  Gulf  ranks 
first  with  63.3%  of  all  AERONET  AOD  measurements  used  in  this  thesis  originating  in 
this  basin.  West  Africa  is  the  source  of  22%  of  AERONET  AOD  measurements, 
followed  by  the  Japan  and  Korea  basins  at  1 1.5%  and  3.2%,  respectively.  Therefore,  the 
composite  climatology  is  heavily  weighted  towards  the  Persian  Gulf  basin.  However, 
even  with  most  of  the  AERONET  AOD  measurements  for  the  2006-2007  period 
originating  in  the  Persian  Gulf  basin,  the  West  Africa  basin  carries,  on  average,  the 
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heaviest  aerosol  loading.  The  Persian  Gulf  basin  ranks  second,  followed  by  the  Korea 
and  Japan  basins.  The  following  sections  will  examine  the  details  behind  the  aerosol 
loading  for  each  basin. 

2.  Persian  Gulf 

This  thesis  includes  AOD  data  from  all  three  Persian  Gulf  AERONET  sites 
(Dhabi,  Dhadnah,  and  Solar  Village)  to  produce  an  AERONET  AOD  climatology  for  the 
Persian  Gulf  basin.  If  this  study  only  includes  data  from  Dhabi  and  Dhadnah,  UAE,  the 
AOD  climatology  will  be  heavily  biased  towards  conditions  at  the  southern  end  of  the 
Persian  Gulf.  Solar  Village,  Saudi  Arabia  is  northwest  of  Dhabi  and  Dhadnah,  and  by 
including  data  from  Solar  Village,  the  AOD  climatology  becomes  more  representative  of 
the  Persian  Gulf  basin  as  a  whole. 

Because  the  AERONET  stations  are  geographically  separated,  there  will  be  some 
subtle  differences  in  AOD  data.  For  example,  while  all  three  AERONET  sites  are  located 
near  dust  source  regions,  only  Dhabi  and  Dhadnah  are  located  near  sources  of  industrial 
pollution  and  maritime  aerosols.  Solar  Village,  on  the  other  hand,  is  located  inland  on  the 
Arabian  Peninsula  (refer  to  Figure  3),  generally  upwind  from  sources  of  industrial 
pollution  (e.g.,  Riyadh)  and  maritime  aerosols.  Consequently,  on  an  annual  basis,  Solar 
Village  should  have  lower  mean  AOD  values  than  Dhabi  and  Dhadnah. 

Figure  13  is  the  two-year  AERONET  0.50pm  AOD  distribution  for  the  Persian 
Gulf  basin  for  2006-2007  using  the  three  stations  Solar  Village,  Dhabi  and  Dhadnah.  The 
figure  contains  both  a  histogram  of  AOD  values  and  cumulative  frequency  distribution. 
Half  of  the  AOD  values  lie  between  0.20  and  0.50,  with  the  median  of  the  distribution  at 
0.35.  In  75%  of  the  cases,  the  aerosol  loading  may  be  light  enough  to  allow  for  scene 
correction,  using  an  AOD  of  0.50  as  a  baseline. 
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Figure  13.  Two-year  distribution  of  AERONET  AOD  (0.50  (am)  values  for  stations 
Solar  Village,  Dhabi  and  Dhadnah  in  the  Persian  Gulf  basin,  for  2006-2007. 

The  histogram  of  AOD  values  (blue  bars)  and  cumulative  frequency 
distribution  (red  curve)  are  denoted.  The  sample  size  is  listed  in  the  upper 
right  corner.  The  vertical  lines  (black),  from  left  to  right,  denote  the  quartiles 
of  the  distribution,  25%  (dashed),  50%  (solid),  75%  (dashed). 

Figure  14  illustrates  the  differences  between  Persian  Gulf  basin  stations  for  2006- 
2007.  The  cumulative  frequency  curve  for  each  individual  station  is  compared  with  the 
cumulative  frequency  curve  for  the  combined  dataset  of  the  three  stations.  On  an  annual 
basis,  aerosol  concentration  tends  to  be  less  at  Solar  Village  than  at  Dhabi  and  Dhadnah. 
This  may  occur  because  Solar  Village  is  located  upwind  of  industrial  pollution  sources 
and  also  sees  less  maritime  aerosols.  There  is  no  significant  difference  in  the  annual 
distribution  of  AERONET  AOD  values  for  Dhabi  and  Dhadnah,  most  likely  due  to 
relatively  close  proximity  to  each  other.  Effects  from  the  Shamal  are  likely  more 
pronounced  over  the  Persian  Gulf  than  inland  over  the  Arabian  Peninsula,  offering 
another  explanation  why  aerosol  concentration  tends  to  be  less  over  Solar  Village. 
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Figure  14.  Two-year  annual  cumulative  frequency  curves  for  the  Persian  Gulf  Basin  for 
2006-2007.  The  cumulative  frequency  curve  for  each  individual  station  is 
compared  with  the  cumulative  frequency  curve  computed  from  the  histogram 
distribution  for  the  combined  dataset  of  the  three  stations  (Persian  Gulf  basin, 

red). 

To  best  illustrate  seasonal  differences  in  AOD  distributions,  this  thesis  divides 
each  year  into  four  seasons,  each  lasting  three  months.  The  seasons  boundaries  are: 
Winter  (December-February),  Spring  (March-May),  Summer  (June-August),  and  Fall 
(September-November).  Figure  15  contains  the  Persian  Gulf  basin  histograms  of 
AERONET  AOD  distributions  by  season  for  2006-2007.  Figure  15a  represents  Winter, 
15b  is  Spring,  15c  is  Summer,  and  15d  is  Fall.  Figures  15a  and  15d  indicate  that  aerosol 
concentrations  are  lowest  over  the  Persian  Gulf  basin  from  September  through  February 
(Fall  and  Winter),  with  approximately  92.5%  of  cases  having  sufficiently  low  aerosol 
loading  to  allow  for  scene  correction.  Figures  15b  and  15c  suggest  that  aerosol 
concentrations  are  highest  from  March  through  August  (Spring  and  Summer),  with 
approximately  68%  of  cases  having  sufficiently  low  aerosol  loading  to  allow  for  scene 
correction  (with  0.50  as  the  baseline  for  correctability).  The  results  from  Figure  15  are 
consistent  with  the  Fourteenth  Weather  Squadron  climatology  for  the  Persian  Gulf  basin 
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(Chapter  II,  Section  B,  Para.  1),  which  found  that  dust  storms  are  more  frequent  in  Spring 
and  Summer  than  in  Winter  and  Fall.  Consequently,  image  analysts  should  expect  more 
uncorrectable  images  in  Spring  and  Summer  versus  Fall  and  Winter. 
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Figure  15.  Seasonal  Distribution  of  AERONET  AOD  (0.50  pm)  values  for  the  Persian 

Gulf  basin  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  are  denoted.  The  sample  size  is 
listed  in  the  upper  right  comer.  Seasons  are  denoted  as  (a)  Winter  (December- 
February),  (b)  Spring  (March-May),  (c)  Summer  (June-August),  and  (d)  Fall 

( S  eptember-N  o  vember) . 

3.  East  Asia 


A  single  basin  AERONET  AOD  climatology  for  East  Asia  is  impractical  for  two 
reasons.  First,  there  is  a  large  geographic  separation  between  western  Korea  and 
southeastern  Japan.  Second,  Korea  is  much  closer  than  Japan  to  dust  and  pollution 
source  regions  in  China  and  Mongolia  (see  Figure  4),  leading  to  the  hypothesis  that,  on 
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an  annual  basis,  AERONET  AOD  measurements  will  be  higher  over  Korea  than  Japan. 
As  a  result,  two  basin  climatologies  are  created,  one  for  Korea,  and  another  for  Japan. 
Only  one  station  (Anmyon)  is  used  in  creating  the  aerosol  climatology  for  Korea,  while 
two  stations  are  used  for  Japan  (Osaka  and  Shirahama). 

Figure  16  is  the  two-year  AERONET  AOD  (0.50pm)  distribution  for  Korea  for 
2006-2007  using  Anmyon.  Figure  17  illustrates  the  same  for  Japan,  using  Shirahama  and 
Osaka.  For  Korea,  fifty  percent  of  the  AOD  values  lie  between  0.14  and  0.49,  with  the 
median  of  the  distribution  at  0.26.  In  75%  of  the  cases  in  Korea,  the  aerosol  loading  may 
be  light  enough  to  allow  for  scene  correction,  using  an  AOD  of  0.50  as  a  baseline. 
Meanwhile,  for  Japan,  50%  of  the  AOD  values  he  between  0.14  and  0.36,  with  the 
median  of  the  distribution  at  0.21.  Additionally,  in  Japan,  the  aerosol  may  be  sufficiently 
light  to  allow  for  scene  correction  in  approximately  87%  of  the  cases.  Furthermore, 
Figures  16  and  17  also  suggest  that,  on  an  annual  basis,  Korea  has  higher  aerosol 
concentrations  than  Japan.  Therefore,  image  analysts  will  have  more  uncorrectable 
imagery  over  Korea  than  Japan. 
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Two-year  distribution  of  AERONET  AOD  (0.50pm)  values  for  Anmyon  in 
the  Korea  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars) 
and  cumulative  frequency  distribution  (red  curve)  are  denoted.  The  sample 
size  is  listed  in  the  upper  right  corner.  The  vertical  lines  (black),  from  left  to 
right,  denote  the  quartiles  of  the  distribution,  25%  (dashed),  50%  (solid),  75% 

(dashed). 


AERONET  AOD  (0.50  (im) 


E 

o 


Figure  17.  Same  as  Figure  16,  but  for  Japan. 
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Figure  18  illustrates  the  differences  between  the  AERONET  AOD  distributions 
for  the  Korea  and  Japan  basins  for  2006-2007.  The  cumulative  frequency  curve  for 
Korea  basin  is  compared  with  the  cumulative  frequency  curve  from  Japan  (see  Figures  16 
and  17).  Each  cumulative  frequency  curve  is  computed  from  the  histogram  distribution 
for  the  each  basin.  Figure  18  indicates  that  the  AERONET  AOD  distributions  are  very 
different  for  Korea  and  Japan,  with  Korea  clearly  carrying  a  heavier  aerosol  load  than 
Japan.  For  example,  for  Korea,  74%  of  the  AERONET  AOD  distribution  has  an  AOD 
less  than  0.50.  Flowever,  for  Japan,  86%  of  the  AERONET  AOD  distribution  has  an 
AOD  less  than  0.50.  Additionally,  the  Korea  basin  has  7%  of  the  AOD  distribution 
greater  than  1.0,  while  the  Japan  basin  only  has  1%  of  the  AOD  distribution  greater  than 
1.0.  Therefore,  the  Korea  and  Japan  AERONET  AOD  datasets  cannot  be  combined  to 
produce  one  single  East  Asia  basin  AERONET  AOD  climatology. 


Figure  18.  Two-year  annual  cumulative  frequency  curves  for  Korea  (blue)  and  Japan 

(red)  for  2006-2007. 

Figure  19  illustrates  the  differences  between  Japan  basin  stations  for  2006-2007. 
There  are  differences  in  the  annual  distribution  of  AERONET  AOD  values  between 
Osaka  and  Shirahama.  Osaka  is  located  next  to  a  sheltered  bay,  surrounded  by  mountains 
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on  three  sides,  while  Shirahama  is  directly  on  the  Pacific  coast.  However,  Shirahama 
may  be  exposed  to  more  dust/pollution  from  sources  in  China,  since  there  are  no 
mountains  that  block  flow  from  China,  and  also  is  likely  sees  more  maritime  aerosols 
than  Osaka.  Therefore,  it  stands  to  reason  that  aerosol  loading  over  Shirahama  is  greater 
than  Osaka.  Also,  many  fewer  AERONET  AOD  measurements  are  available  for  Osaka 
than  Shirahama,  which  explains  why  the  Japan  basin  composite  cumulative  frequency 
curve  closely  mirrors  the  Shirahama  curve. 


Figure  19.  Same  as  Figure  18,  but  for  the  Japan  basin. 

Figure  20  contains  the  seasonal  distributions  of  AERONET  AOD  values  for 
Korea.  Figures  20a  and  20d  indicate  that  aerosol  concentrations  are  lowest  over  Korea 
from  September  through  February  (Fall  and  Winter),  with  approximately  90%  of  cases 
having  sufficiently  low  aerosol  loading  to  allow  for  scene  correction.  Figures  20b  and 
20c  suggest  that  aerosol  concentrations  are  highest  from  March  through  August  (Spring 
and  Summer),  with  approximately  58%  of  cases  having  sufficiently  low  aerosol  loading 
to  allow  for  scene  correction.  The  results  from  Figure  20  agree  with  the  Fourteenth 
Weather  Squadron  climatology  (Chapter  II,  Section  B,  Para.  2),  that  dust  storms  are  more 
frequent  in  Spring  and  Summer  than  in  Winter  and  Fall. 
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Figure  20.  Seasonal  Distribution  of  AERONET  AOD  (0.50  pm)  values  for  the  Korea 
basin  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  are  denoted.  The  sample  size  is 
listed  in  the  upper  right  comer.  Seasons  are  denoted  as  (a)  Winter  (December- 
February),  (b)  Spring  (March-May),  (c)  Summer  (June-August),  and  (d)  Fall 

( S  eptember-N  o  vember) . 

Figure  21  contains  the  seasonal  distributions  of  AERONET  AOD  values  for 
Japan.  Figures  21a  and  2 Id  indicate  that  aerosol  concentrations  are  lowest  over  Japan 
from  September  through  February  (Fall  and  Winter),  with  more  than  95%  of  cases  having 
sufficiently  low  aerosol  loading  to  allow  for  scene  correction  (with  0.50  as  the  baseline 
for  correctability).  Figures  21b  and  21c  suggest  that  aerosol  concentrations  are  highest 
from  March  through  August  (Spring  and  Summer),  with  approximately  80%  percent  of 
cases  having  sufficiently  low  aerosol  loading  to  allow  for  scene  correction.  The  results 
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from  Figure  21  are  consistent  with  the  Fourteenth  Weather  Squadron  climatology 
(Chapter  II,  Section  B,  Para.  2),  that  dust  storms  are  more  frequent  in  Spring  and  Summer 
than  in  Winter  and  Fall  over  Japan. 
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Figure  2 1 .  Same  as  Figure  20  but  for  the  Japan  basin. 

Additionally,  aerosol  loading  appears  to  be  higher  over  Korea  than  Japan.  As 
Korea  is  closer  to  dust  and  pollution  sources  in  China  and  Mongolia,  Korea  should  have  a 
heavier  aerosol  load  than  Japan.  When  analyzing  imagery  over  East  Asia,  an  analyst 
should  expect  more  aerosol  loading  over  Korean  imagery  versus  Japanese  imagery,  with 
most  uncorrectable  images  occurring  during  Spring  and  Summer. 

4.  West  Africa 

AOD  data  from  both  West  Africa  AERONET  sites  is  used  to  produce  an  AOD 
climatology  for  the  West  African  basin.  If  data  from  only  Ilorin,  Nigeria  are  included,  for 
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example,  the  AOD  climatology  will  be  heavily  biased  towards  conditions  at  the  Ilorin. 
Djougou,  Benin  is  northwest  of  Ilorin,  and  by  including  data  from  Djougou,  the  AOD 
climatology  becomes  more  representative  of  the  West  African  basin  as  a  whole,  which 
helps  verify  the  significant  basin  differences  compared  to  the  Persian  Gulf  and  East  Asia. 

Because  the  AERONET  stations  are  geographically  separated,  there  will  be  some 
subtle  differences  in  AOD  data.  Additionally,  no  AOD  data  is  available  at  0.50pm  for 
Djougou.  Instead,  this  thesis  uses  AERONET  AOD  data  at  0.44pm  for  Djougou.  Since 
AOD  generally  increases  with  decreasing  wavelength,  we  expect  the  Djougou  results  to 
be  biased  toward  slightly  higher  values. 

Figures  22  and  23  are  the  two-year  AERONET  AOD  distributions  for  the  two 
West  Africa  stations.  Overall  aerosol  loading  in  West  Africa  is  quite  high.  For  Djougou, 
half  of  the  AOD  values  lie  between  0.42  and  0.89,  with  the  median  of  the  distribution  at 
0.61.  For  Ilorin,  half  of  the  AOD  values  lie  between  0.40  and  1.05,  with  the  median  of 
the  distribution  at  0.70.  The  aerosol  loading  may  be  light  enough  to  allow  for  scene 
correction  in  only  35%  of  the  cases  in  West  Africa,  using  an  AOD  of  0.50  as  a  baseline. 

Additionally,  Figures  22  and  23  suggest  that  annual  aerosol  concentrations  over 
West  Africa  are  much  higher  than  for  the  Persian  Gulf  and  East  Asia.  For  example,  for 
West  Africa,  22.5%  of  AERONET  AOD  measurements  taken  during  2006-2007  are 
greater  than  1.0,  versus  3%  of  AERONET  AOD  measurements  for  the  Persian  Gulf 
basin.  Aerosol  loading  over  West  Africa  is  greater  than  aerosol  loading  for  the  Persian 
Gulf  and  East  Asia  because  of  the  proximity  of  Djougou  and  Ilorin  to  the  Sahara  Desert. 
Because  both  stations  are  located  in  the  Sahel  region  of  Africa,  Djougou  and  Ilorin  are 
subject  to  strong  Hannattan  winds  that  blow  over  West  Africa  during  the  dry  season 
(roughly  November  through  March).  These  winds  blanket  the  region  with  high 
concentrations  of  dust  and  sand  from  the  Sahara  Desert,  and  are  a  major  source  of 
aerosols  over  West  Africa  (Westphal  1988).  The  other  major  source  of  aerosols  during 
the  dry  season  is  biomass  burning,  generally  advected  from  the  regions  south  and  east  of 
the  West  Africa  stations. 
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Figure  22.  Two-year  distribution  of  AERONET  AOD  (0.44|tim)  values  for  station 
Djougou  in  the  West  Africa  basin,  for  2006-2007.  The  histogram  of  AOD 
values  (blue  bars)  and  cumulative  frequency  distribution  (red  curve)  are 
denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical  lines 
(black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  24  illustrates  the  differences  between  West  Africa  basin  stations  for  2006- 
2007.  At  AODs  below  0.55,  there  is  little  difference  in  AOD  distributions  between 
Djougou  and  Ilorin.  This  implies  that  aerosol  source  differences  are  larger  than  the  effect 
of  the  shorter  wavelength  used  for  Djougou.  At  AODs  above  0.55,  there  is  a  greater 
difference  in  AOT  distributions.  Ilorin  is  much  closer  to  Lake  Chad,  which  is  a  major 
dust  source  region  in  West  Africa.  It  is  possible  that  Ilorin,  given  its  proximity  to  Lake 
Chad,  receives  more  dust  than  Djougou.  Thus,  the  AOD  distribution  for  Ilorin  shows  a 
heavier  aerosol  concentration  as  compared  to  Djougou. 


Figure  24.  Two-year  annual  cumulative  frequency  curves  for  the  West  Africa  basin  for 
2006-2007.  The  Ilorin  (0.50pm,  blue)  and  Djougou  (0.44pm,  green)  curves 

are  denoted. 

Figure  25  contains  seasonal  West  Africa  basin  histograms  for  Ilorin,  which 
illustrate  AERONET  AOD  distributions  by  season  for  2006-2007.  Figures  25c  and  25d 
indicate  that  aerosol  concentrations  are  lowest  over  Ilorin  from  June  through  November 
(wet  season),  with  approximately  75%  of  cases  having  sufficiently  low  aerosol  loading  to 
allow  for  scene  correction,  using  an  AOD  value  of  0.50  as  a  baseline  for  image 
correctability.  Figures  25a  and  25b  suggest  that  aerosol  concentrations  are  highest  from 
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December  through  May  (dry  season),  with  only  approximately  20%  of  cases  having 
sufficiently  low  aerosol  loading  to  allow  for  scene  correction.  The  results  from  Figure  25 
are  consistent  with  the  climatology  produced  by  the  Fourteenth  Weather  Squadron,  that 
dust  storms  are  more  frequent  in  the  dry  season  than  in  the  wet  season  for  the  West 
Africa  basin.  Seasonal  AERONET  AOD  histograms  for  Djougou  are  included  in 
Appendix  B  (Figures  81-84),  and  display  the  same  trends  as  Ilorin. 
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Figure  25.  Seasonal  Distribution  of  AERONET  AOD  (0.50  pm)  values  for  Ilorin  in  the 
West  Africa  basin  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars) 
and  cumulative  frequency  distribution  (red  curve)  are  denoted.  The  sample 
size  is  listed  in  the  upper  right  corner.  Seasons  are  denoted  as  (a)  Winter 
(December-February),  (b)  Spring  (March-May),  (c)  Summer  (June -August), 
and  (d)  Fall  (September-November). 

As  shown  in  Figures  25a  and  25b,  AERONET  AOD  distributions  are  shifted 
towards  the  high  end  of  the  AOD  spectrum  during  the  dry  season  (Winter  and  Spring). 
Based  on  climatology,  West  Africa  stations  should  see  higher  aerosol  concentrations 
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during  the  dry  season,  and  that  is  exactly  what  is  shown  by  the  AERONET  AOD 
distributions.  With  the  approach  of  the  rainy  season  (Summer  and  Fall),  the  Harmattan 
winds  die  down,  and  the  AERONET  AOD  distributions  shift  towards  lower  values.  Also, 
the  number  of  cloudy  days  increases,  reducing  the  frequency  that  AERONET 
measurements  are  taken,  and  thus,  lowering  the  sample  sizes. 

In  summary,  an  analyst  viewing  targets  over  West  Africa  will  have  to  apply  the 
greatest  amount  of  corrections  in  the  dry  season,  due  to  the  increased  frequency  of  dust 
storms  over  the  region.  However,  during  the  rainy  season,  while  there  will  be  more 
correctable  imagery,  there  will  be  fewer  images  to  work  with  because  of  increased  cloud 
cover. 

B.  MODIS/ AERONET  COMPARISON 

1.  Persian  Gulf 

The  purpose  of  the  MODIS/ AERONET  comparison  is  to  establish  that  MODIS- 
retrieved  AOD  values  at  0.55pm  are  representative  of  ground-truth  AOD  values 
throughout  each  geographic  basin.  AOD  data  from  only  two  Persian  Gulf  AERONET 
sites  are  included  in  the  MODIS/AERONET  comparison  for  the  Persian  Gulf  basin. 
MODIS -retrieved  AOD  values  are  from  over  the  ocean  only.  Solar  Village,  Saudi  Arabia 
is  approximately  417  km  inland  from  the  Persian  Gulf,  so  there  are  no  MODIS -retrieved 
AOD  values  in  the  vicinity  of  Solar  Village.  Consequently,  AERONET  AOD  data  from 
Solar  Village  are  not  used  in  the  MODIS/AERONET  comparison.  The  Dhabi,  UAE  and 
Dhadnah,  UAE  AERONET  sites  are  both  located  along  the  coast,  so  this  study  is  able  to 
make  point-to-point  MODIS/AERONET  AOD  comparisons  for  both  sites  (see  Figure  6). 
The  datasets  from  Dhabi  and  Dhadnah  are  combined  to  produce  the  Persian  Gulf  basin 
MODIS/AERONET  AOD  comparison. 

MODIS  and  AERONET  AOD  values  need  to  be  coincident  in  time  and  space  for 
proper  comparison.  The  MODIS  grid  point  closest  to  each  AERONET  site  is  used  (see 
Figure  6).  Each  MODIS  AOD  value  is  a  daily  mean  of  the  AOD  values  for  a  1°  x  1°  grid 
box.  MODIS  overpass  times  for  the  Persian  Gulf  basin  are  approximately  1030LT 
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(0630UTC)  and  1330LT  (0930UTC).  For  each  AERONET  site,  AOD  data  is  averaged 
between  1000LT  (0600UTC)  and  1400LT  (1000UTC)  to  be  consistent  with  the  time  span 
of  the  MODIS  satellite  passes.  On  each  day  when  both  MODIS  and  AERONET  AOD 
values  are  available  for  the  time  window,  a  one-to-one  comparison  of  AOD  values  is 
performed. 

Figure  26  contains  a  scatter  plot  of  AERONET  and  MODIS  AOD  values  for  the 
Persian  Gulf  basin  for  the  two-year  period  from  2006  through  2007.  Agreement  between 
the  MODIS  and  AERONET  values  is  quite  good.  The  correlation  coefficient  is  0.90,  and 
the  scatter  does  not  imply  bias  in  the  comparison.  At  an  AOD  greater  than  0.50,  the 
spread  increases  slightly  with  increasing  AOD.  The  spread  may  increase  because  of 
greater  differences  in  time  and  space  between  the  observation  techniques. 


Figure  26.  Annual  scatter  plot  of  daily  AERONET  and  MODIS  AOD  values  for  the 
Persian  Gulf  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one- 
to-one  line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size 
(N)  are  listed  in  the  upper  left  comer. 

Figure  27  contains  scatter  plots  of  AERONET  and  MODIS  AOD  values  by 

season  for  the  2006-2007  timeframe.  There  is  good  correlation  for  all  seasons.  As 

before,  the  spread  increases  as  AOD  increases.  Spring  and  Summer  display  higher 
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average  AOD  values,  reinforcing  that  Spring  and  Summer  are  the  seasons  with  higher 
aerosol  loading  over  the  Persian  Gulf  basin.  The  implications  are  that  for  Spring  and 
Summer  over  the  Persian  Gulf  basin,  an  analyst  should  have  a  higher  percentage  of 
images  that  are  uncorrectable.  Additionally,  Fall  and  Winter  are  the  seasons  of  the  year 
that  may  allow  for  more  images  to  be  corrected.  Finally,  Figure  27  confirms  that,  on  a 
seasonal  basis,  MODIS-retrieved  AOD  values  are  representative  of  the  land-based 
AERONET  AOD  values  for  the  Persian  Gulf  basin. 


Figure  27.  Persian  Gulf  Basin  AERONET  AOD  (0.50pm)  versus  MODIS  AOD 

(0.55pm)  scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c) 
Summer,  (d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line 
(blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are 

listed  in  the  upper  left  comer. 

Figure  28  displays  MODIS  and  AERONET  annual  cumulative  frequency  curves 
from  the  Persian  Gulf  basin  for  the  two-year  period  (2006-2007).  The  MODIS 
cumulative  frequency  curve  is  derived  from  the  two-year  Persian  Gulf  basin  MODIS 
AOD  distribution  (0.55pm)  histogram,  located  in  Appendix  B  (Figure  85).  The 
AERONET  cumulative  frequency  curve  is  taken  from  the  two-year  basin  AERONET 
AOD  distribution  (0.50pm)  histogram  (Figure  13).  The  MODIS  cumulative  frequency 
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curve  closely  mirrors  the  AERONET  AOD  curve,  suggesting  that  on  an  annual  basis 
there  is  little  difference  in  the  distribution  of  AOD  values  for  both  MODIS  and 
AERONET.  For  example,  40%  of  both  MODIS  and  AERONET  AOD  values  are  less 
than  an  AOD  of  0.30.  The  region  of  most  spread  between  the  MODIS  and  AERONET 
curves  is  between  AODs  of  0.05  to  0.25,  with  the  MODIS  curve  indicating  slightly 
higher  AODs  than  AERONET.  At  AODs  from  0.05  to  0.25,  MODIS  portrays  heavier 
aerosol  loading  compared  to  what  is  actually  occurring.  Turbid  waters  in  the  Persian 
Gulf  add  TOA  radiance,  influencing  the  MODIS  AOD  retrieval,  resulting  in  an 
overestimate  of  the  AOD  values  over  water.  For  optical  depths  above  about  0.25,  aerosol 
radiance  contribution  may  dominate  the  contribution  from  the  surface. 


Figure  28.  Annual  cumulative  frequency  curves  for  the  Persian  Gulf  basin  MODIS  (red) 

and  AERONET  (blue)  AOD  distributions. 

Table  1  provides  a  summary  of  the  statistical  results  for  the  MODIS/AERONET 
AOD  comparison  for  the  Persian  Gulf  basin.  Sample  sizes  are  sufficiently  large  such  that 
valid  conclusions  can  be  drawn  from  the  data.  Correlation  coefficients  range  from  0.87 
to  0.90,  which  suggest  a  high  degree  of  correlation  between  MODIS  and  AERONET 
AOD  values.  Thus,  the  implication  is  that  the  MODIS  AOD  values  for  the  Persian  Gulf 
basin  can  be  used  with  a  high  degree  of  confidence  in  areas  lacking  AERONET  sites. 


40 


Table  1.  MODIS/AERONET  AOD  Correlations  for  Persian  Gulf  basin  for  2006-2007 


N  (sample  size) 

0.90 

743 

0.88 

140 

Spring 

0.87 

218 

Summer 

0.87 

248 

Fall 

0.87 

137 

In  summary,  the  MODIS  AOD  retrieval  at  0.55pm  can  provide  reliable  estimates 
of  AOD  values  in  the  Persian  Gulf  basin  for  areas  not  located  near  AERONET  sites. 
While  the  AERONET  sites  used  here  are  along  coasts,  they  are  still  located  on  land  and 
are  located  closer  to  dust  and  pollution  sources.  The  proximity  to  the  aerosol  sources 
may  account  for  some  of  the  slight  differences  between  MODIS  and  AERONET  AOD 
measurements  for  the  basin.  Additionally,  local  weather  conditions  at  the  AERONET 
sites  (i.e.,  sea  and  land  breezes)  no  doubt  affect  aerosol  loading.  In  turn,  some  of  the 
differences  between  MODIS  and  AERONET  AOD  measurements  may  come  about  due 
to  local  weather  conditions  at  the  AERONET  sites. 

2.  East  Asia 

a.  Korea 

Only  data  from  the  Anmyon,  South  Korea  AERONET  site  is  included  in 
the  MODIS/AERONET  comparison  for  the  Korea  basin.  AOD  values  from  the  MODIS 
grid  point  closest  to  Anmyon  is  used  in  this  comparison  (see  Figure  7).  As  previously 
described,  AERONET  AOD  values  are  averaged  over  a  time  span  to  match  the  time 
frame  of  the  mean  MODIS  AOD  value.  The  MODIS  satellites  image  the  Korea  basin  at 
approximately  1030LT  (0130UTC)  and  1330LT  (0430UTC).  The  mean  AERONET 
AOD  value  is  computed  for  the  time  period  from  1000-1400LT  (0100-0500UTC),  for 
each  cloud-free  day  for  which  MODIS  AOD  values  are  available. 
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Figure  29  contains  a  scatter  plot  of  AERONET  and  MODIS  AOD  values 
for  the  Korea  basin  for  the  two-year  period  from  2006  through  2007.  Agreement 
between  the  MODIS  and  AERONET  values  is  quite  good.  The  scatter  does  not  imply 
bias  in  the  comparison,  and  the  correlation  coefficient  is  0.91. 


Figure  29.  Annual  scatter  plot  of  daily  AERONET  and  MODIS  AOD  values  for  the 
Korea  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one-to-one 
line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are 

listed  in  the  upper  left  corner. 

Figure  30  contains  scatter  plots  of  AERONET  and  MODIS  AOD  values 
by  season  for  the  2006-2007  time  period.  The  correlation  is  high  for  all  seasons.  Spring 
and  Summer  display  higher  average  AOD  values,  reinforcing  that  Spring  and  Summer 
are  the  seasons  with  higher  aerosol  loading  over  the  Korea  basin.  The  implications  are 
that  for  Spring  and  Summer  over  the  Korea  basin,  an  analyst  should  have  a  higher 
percentage  of  images  that  are  uncorrectable,  with  fewer  correctable  images. 
Additionally,  Fall  and  Winter  are  the  seasons  of  the  year  that  may  allow  for  more  images 
to  be  corrected. 
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Figure  30  also  shows  that  sample  size  is  small  for  the  Korea  basin.  The 
small  sample  size  serves  to  increase  the  uncertainty  behind  the  MODIS  AOD  retrieval. 
One  possible  reason  for  small  sample  size  is  the  amount  of  cloud  cover  over  Korea.  In 
comparison  to  the  Persian  Gulf,  Korea  is,  on  average,  much  cloudier,  resulting  in  fewer 
days  on  which  MODIS  and  AERONET  AOD  retrievals  can  be  made.  Equipment  shutoff 
or  failure  may  also  contribute  to  the  small  sample  size  for  the  Korea  basin. 


Figure  30.  Korea  Basin  AERONET  AOD  (0.50pm)  versus  MODIS  AOD  (0.55pm) 
scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c)  Summer, 
(d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line  (blue)  are 
denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are  listed  in  the 

upper  left  corner. 

Figure  3 1  displays  MODIS  and  AERONET  annual  cumulative  frequency 
curves  from  the  Korea  basin  for  the  two-year  period  (2006-2007).  The  MODIS 
cumulative  frequency  curve  is  derived  from  the  two-year  Korea  basin  MODIS  AOD 
distribution  (0.55pm)  histogram,  located  in  Appendix  B  (Figure  87).  The  AERONET 
cumulative  frequency  curve  is  taken  from  the  two-year  basin  AERONET  AOD 
distribution  (0.50pm)  histogram  (Figure  16).  The  MODIS  curve  is  far  to  the  right  of  the 
AERONET  curve,  indicating  that  on  an  annual  basis,  MODIS  displays  a  tendency  to 
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portray  heavier  aerosol  loading  compared  to  AERONET.  Turbid,  shallow  waters  in  the 
Yellow  Sea  and  a  small  AERONET  dataset  for  Korea  may  explain  why  MODIS  displays 
a  tendency  to  overestimate  aerosol  loading. 


Figure  3 1 .  Annual  cumulative  frequency  curves  for  the  Korea  basin  MODIS  (red)  and 

AERONET  (blue)  AOD  distributions. 

Table  2  provides  a  summary  of  the  statistical  results  for  the 
MODIS/AERONET  AOD  comparison  for  the  Korea  basin.  Correlation  coefficients 
range  from  0.91  to  0.94.  Table  2,  in  addition  to  Figures  29-31;  indicate  that  the  MODIS 
AOD  retrieval  at  0.55pm  may  provide  close  to  ground  truth  AOD  values  in  the  Korea 
basin  for  areas  not  located  near  AERONET  sites.  Larger  AERONET  and  MODIS  AOD 
sample  sizes  will  help  narrow  down  the  uncertainty  as  to  whether  or  not  MODIS  can 
reliably  provide  AOD  measurements  over  the  Korea  basin. 
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Table  2.  MODIS/AERONET  AOD  Correlations  for  Korea  basin  for  2006-2007 


R  (correlation) 

N  (sample  size) 

Annual 

0.91 

38 

Winter 

0.94 

15 

Spring 

0.94 

7 

Summer 

0.92 

7 

Fall 

0.88 

9 

b.  Japan 

Shirahama  MODIS  AOD  data  are  included  in  the  MODIS/AERONET 
comparison  for  the  Japan  basin.  Osaka  MODIS  AOD  data  are  not  included  because 
limited  MODIS  data  are  available  in  the  proximity  of  Osaka.  AOD  values  from  the 
MODIS  grid  point  closest  to  Shirahama  are  used  in  this  comparison  (see  Figure  7).  As 
previously  described,  AERONET  AOD  values  are  averaged  over  time  to  match  the  mean 
MODIS  AOD  values.  The  MODIS  satellites  image  the  Japan  basin  at  approximately 
1030L  (0130UTC)  and  1330L  (0430UTC).  The  mean  AERONET  AOD  value  is 
computed  for  1000-1400L  (0100-0500UTC),  for  each  cloud-free  day  for  which  MODIS 
AOD  values  are  available. 

Figure  32  contains  a  scatter  plot  of  AERONET  and  MODIS  AOD  values 
for  the  Japan  basin  for  the  two-year  period  (2006-2007).  Agreement  between  the  MODIS 
and  AERONET  values  is  very  good,  with  the  correlation  coefficient  equal  to  0.86. 
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Figure  32.  Annual  scatter  plot  of  daily  AERONET  and  MODIS  AOD  values  for 
Shirahama  in  the  Japan  basin  for  2006-2007.  The  daily  AOD  values  (red 
stars)  and  one-to-one  line  (blue)  are  denoted.  The  correlation  coefficient  (R) 
and  sample  size  (N)  are  listed  in  the  upper  left  corner. 

Figure  33  contains  scatter  plots  of  AERONET  and  MODIS  AOD  values 
by  season  for  2006-2007.  The  correlation  is  not  consistent  across  the  seasons.  One 
possible  reason  may  be  the  low  sample  size.  Further  explanation  of  the  differences  in 
correlation  by  season  is  beyond  the  scope  of  this  thesis. 
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Figure  33.  Japan  Basin  AERONET  AOD  (0.50pm)  versus  MODIS  AOD  (0.55pm) 
scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c)  Summer, 
(d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line  (blue)  are 
denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are  listed  in  the 

upper  left  corner. 

Figure  34  displays  MODIS  and  AERONET  annual  cumulative  frequency 
curves  from  the  Japan  basin  for  the  two-year  period  (2006-2007).  The  MODIS 
cumulative  frequency  curve  is  derived  from  the  two-year  Japan  basin  MODIS  AOD 
distribution  (0.55pm)  histogram,  located  in  Appendix  B  (Figure  89).  The  AERONET 
cumulative  frequency  curve  is  taken  from  the  two-year  basin  AERONET  AOD 
distribution  (0.50pm)  histogram  (Figure  17).  The  MODIS  cumulative  frequency  curve, 
in  general,  closely  mirrors  the  AERONET  AOD  curve,  suggesting  that  on  an  annual  basis 
there  are  only  slight  differences  in  the  distribution  of  AOD  values  for  both  MODIS  and 
AERONET.  Below  an  AOD  of  0.30,  the  MODIS  curve  is  slightly  to  the  right  of  the 
AERONET  curve,  implying  that  MODIS  slightly  overestimates  AOD  values  when  the 
actual  AOD  is  less  than  or  equal  to  0.30  (possibly  because  of  turbid  waters).  Above  an 
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AOD  of  0.30,  the  MODIS  curve  is  slightly  to  the  left  of  the  AERONET  curve,  suggesting 
that  MODIS  slightly  underestimates  AOD  values  when  the  actual  AOD  is  less  than  or 
equal  to  0.30. 


Figure  34.  Annual  cumulative  frequency  curves  for  the  Japan  basin  MODIS  (red)  and 

AERONET  (blue)  AOD  distributions. 

Table  3  summarizes  the  statistical  results  for  the  MODIS/AERONET 
AOD  comparison  for  the  Japan  basin.  Correlation  coefficients  range  from  0.76  to  0.94, 
which  suggest  a  high  degree  of  correlation  between  MODIS  and  AERONET  AOD 
values. 

Table  3.  MODIS/AERONET  AOD  Correlations  for  Japan  basin  for  2006-2007 


R  (correlation) 

N  (sample  size) 

Annual 

0.86 

237 

Winter 

0.94 

73 

Spring 

0.76 

67 

Summer 

0.89 

58 

Fall 

0.78 

39 
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In  summary,  like  the  Persian  Gulf  basin,  the  MODIS  AOD  retrieval  at 
0.55  pm  can  provide  a  good  estimate  of  AOD  values  in  the  Japan  basin  for  areas  not 
located  near  AERONET  sites.  While  the  AERONET  site  used  in  this  thesis  is  along  the 
coast,  the  site  is  still  located  on  land,  and  is  located  closer  to  dust  and  pollution  sources. 
The  proximity  to  the  aerosol  sources  may  account  for  some  of  the  slight  differences 
between  MODIS  and  AERONET  AOD  measurements  for  the  basin.  Additionally,  local 
weather  conditions  at  the  AERONET  site  (i.e.,  sea  and  land  breezes)  no  doubt  affect 
aerosol  loading.  In  turn,  some  of  the  differences  between  MODIS  and  AERONET  AOD 
measurements  may  come  about  due  to  local  weather  conditions  at  the  AERONET  sites. 

3.  West  Africa 

Only  AOD  data  from  Ilorin,  Nigeria  are  included  in  the  MODIS/AERONET 
comparison  for  the  West  Africa  basin.  The  MODIS  AOD  retrieval  is  an  over  water 
retrieval.  A  more  accurate  MODIS/AERONET  comparison  will  result  using  only  Ilorin 
because  Djougou,  Benin  is  more  than  167  km  further  inland  than  Ilorin.  AOD  values 
from  the  MODIS  grid  point  closest  to  Ilorin  are  used  in  this  comparison  (see  Figure  8). 
As  previously  described,  AERONET  AOD  values  are  averaged  over  time  to  match  the 
mean  MODIS  AOD  values.  The  MODIS  satellites  image  the  West  Africa  basin  at 
approximately  1030L  (0930UTC)  and  1330L  (1230UTC).  The  mean  AERONET  AOD 
value  is  computed  for  1000-1400L  (0900-1300UTC),  for  each  cloud-free  day  for  which 
MODIS  AOD  values  are  available. 

Figure  35  contains  a  scatter  plot  of  AERONET  and  MODIS  AOD  values  for  the 
West  Africa  basin  for  the  two-year  period  from  2006  through  2007.  Agreement  between 
the  MODIS  and  AERONET  values  is  not  nearly  as  good  for  West  Africa  as  it  is  for  the 
Persian  Gulf  and  Japan  basins.  MODIS  AOD  values  have  a  significant  low  bias.  The 
source  of  the  low  bias  may  be  due  to  the  distance  between  the  Ilorin  AERONET  station 
and  the  coast.  Ilorin  is  over  217  km  inland,  and  is  much  closer  to  dust  and  pollution 
sources  than  the  over-ocean  values  of  the  MODIS  AOD  retrieval.  Consequently,  AOD 
values  measured  over  Ilorin  should  be  greater  than  over-ocean  MODIS  AOD  values. 
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Figure  35.  Annual  scatter  plot  of  daily  AERONET  and  MODIS  AOD  values  for  Ilorin 
in  the  West  Africa  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and 
one-to-one  line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample 
size  (N)  are  listed  in  the  upper  left  corner. 

Figure  36  contains  scatter  plots  of  AERONET  and  MODIS  AOD  values  by 
season  for  the  2006-2007  timeframe.  A  low  bias  exists  in  the  seasonal  comparisons.  The 
MODIS  retrieval  regularly  underestimates  AOD  values,  especially  in  the  dry  season 
(Winter  and  Spring).  Sample  size  is  small  during  Summer  and  Fall.  The  Summer  and 
Fall  encompasses  the  cloudy,  wet  season  in  West  Africa,  so  the  number  of  cloud-free 
days  on  which  MODIS  can  make  retrievals  is  low.  Even  with  a  low  sample  size,  Figure 
36  suggests  that  the  MODIS  AOD  retrieval  may  not  provide  representative  AOD  values 
for  the  West  Africa  basin,  due  to  the  proximity  of  the  AERONET  sites  to  the  dust  and 
pollution  sources. 
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Figure  36.  West  Africa  Basin  AERONET  AOD  (0.50pm)  versus  MODIS  AOD 

(0.55 |iun)  scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c) 
Summer,  (d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line 
(blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are 

listed  in  the  upper  left  comer. 

Figure  37  displays  MODIS  and  AERONET  annual  cumulative  frequency  curves 
from  the  West  Africa  basin  for  the  two-year  period  (2006-2007).  The  MODIS 
cumulative  frequency  curve  is  derived  from  the  two-year  West  Africa  basin  MODIS 
AOD  distribution  histogram,  located  in  Appendix  B  (Figure  91).  The  AERONET 
cumulative  frequency  curve  is  taken  from  the  two-year  basin  AERONET  AOD 
distribution  histogram  (Figure  23).  The  MODIS  cumulative  frequency  curve  remains  to 
the  left  of  the  AERONET  curve,  suggesting  that  on  an  annual  basis  the  MODIS  AOD 
retrieval  underestimates  aerosol  loading  over  West  Africa.  Due  to  the  distance  between 
the  MODIS  grid  point  and  aerosol  sources,  it  is  reasonable  that  on  an  annual  basis,  the 
MODIS  AOD  retrieval  underestimates  aerosol  loading.  The  implication  from  Figure  37 
is  that  analysts  should  expect  MODIS  to  underestimate  aerosol  loading  over  the  West 
Africa  basin. 
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Figure  37.  Annual  cumulative  frequency  curves  for  the  West  Africa  basin  MODIS  (red) 
and  AERONET  (blue)  AOD  distributions.  The  y-axis  indicates  the 
cumulative  probability  of  occurrence  of  an  AOD.  The  x-axis  indicates  the 

AOD  values. 


Table  4  provides  a  summary  of  the  statistical  results  for  the  MODIS/AERONET 
AOD  comparison  for  the  West  Africa  basin.  Correlation  coefficients  range  from  0.69  to 
0.92.  MODIS  AOD  results  from  the  Fall  season  appear  to  be  more  reliable  than  for  the 
other  seasons. 

Table  4.  MODIS/AERONET  AOD  Correlations  for  West  Africa  basin  for  2006-2007 


R  (correlation) 

N  (sample  size) 

Annual 

0.77 

107 

Winter 

0.69 

71 

Spring 

0.80 

24 

Summer 

N/A 

1 

Fall 

0.92 

11 

In  summary,  MODIS-retrieved  AOD  values  consistently  underestimate  aerosol 
loading  throughout  the  West  Africa  basin.  As  such,  over-water  MODIS-retrieved  AOD 
values  are  not  a  reliable  proxy  for  AERONET  AOD  values  over  West  Africa. 
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4.  Summary 


The  MODIS  AOD  retrieval  at  0.5 5 jam  consistently  provides  accurate  AOD  values 
for  the  Persian  Gulf  and  Japan  basins.  The  MODIS  AOD  retrieval  is  not  as  accurate  for 
the  Korea  and  West  Africa  basins.  The  disparity  between  MODIS-retrieved  and 
AERONET  AOD  values  over  Korea  may  be  due  to  the  small  sample  sizes  over  the 
Korean  Peninsula.  West  Africa,  on  the  other  hand,  has  more  MODIS  and  AERONET 
AOD  data.  However,  the  MODIS  results  are  verified  against  an  AERONET 
approximately  217  km  from  the  Atlantic  coast,  which  naturally  introduces  error  into  the 
comparison. 

C.  NAAPS/AERONET  COMPARISON 

1.  Persian  Gulf 

The  purpose  of  the  NAAPS/AERONET  comparison  is  to  consider  if  the  NAAPS- 
modeled  AOD  values  at  0.55pm  are  representative  of  ground-based  AOD  measurements 
throughout  each  geographic  basin.  Because  NAAPS  AOD  values  are  provided  only  over 
water,  no  AERONET  data  from  Solar  Village,  Saudi  Arabia  are  included.  The  NAAPS 
AOD  grid  point  closest  to  each  coastal  station,  Dhabi  and  Dhadnah,  UAE,  are  used  for 
the  comparison.  Both  stations  are  combined  to  produce  the  Persian  Gulf  basin 
NAAPS/AERONET  comparison. 

Figure  38  is  the  two-year  NAAPS  0.55pm  AOD  distribution  for  the  Persian  Gulf 
basin  for  2006-2007.  Half  of  the  AOD  values  lie  between  0.25  and  0.45,  with  the  median 
of  the  distribution  at  0.32.  In  nearly  80%  of  the  cases  the  aerosol  loading  may  be  light 
enough  to  allow  for  scene  correction,  using  an  AOD  of  0.50  as  a  baseline.  In 
comparison,  Persian  Gulf  basin  AERONET  AOD  data  (Figure  13)  shows  that  (a)  the 
median  of  the  distribution  is  0.35  and  (b)  75%  of  the  cases  show  light  enough  aerosol 
loading  to  allow  for  scene  correction.  Therefore,  the  comparison  between  Figures  13  and 
38  suggests  that  for  the  two-year  period  (2006-2007),  the  NAAPS  aerosol  loading  over 
the  Persian  Gulf  basin  closely  matches  the  AERONET  aerosol  loading. 
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Figure  38.  Two-year  distribution  of  daily  mean  0.5 5  jam  NAAPS  AOD  values  for 
stations  Dhabi  and  Dhadnah  in  the  Persian  Gulf  basin,  for  2006-2007.  The 
histogram  of  AOD  values  (blue  bars)  and  cumulative  frequency  distribution 
(red  curve)  corresponding  to  the  histogram  are  denoted.  The  sample  size  is 
listed  in  the  upper  right  comer.  The  vertical  lines  (black),  from  left  to  right, 
denote  the  quartiles  of  the  distribution,  25%  (dashed),  50%  (solid),  75% 

(dashed). 

Figure  39  contains  a  scatter  plot  of  AERONET  and  NAAPS  AOD  values  for  the 
Persian  Gulf  basin  for  the  two-year  period  from  2006  through  2007.  Agreement  between 
the  NAAPS  and  AERONET  values  is  good.  The  scatter  implies  no  bias  in  the 
comparison.  The  amount  of  scatter  does  increase  at  higher  AOD  values  (i.e.  above  0.60), 
suggesting  that  the  uncertainty  of  NAAPS  estimates  for  the  Persian  Gulf  increases  with 
AOD.  The  uncertainty  of  NAAPS  estimates  may  increase  with  AOD  due  to  errors  in  the 
modeled  wind  field,  aerosol  type,  and  dust  source.  Additionally,  the  uncertainty  of  the 
MODIS  AOD  retrieval  assimilated  into  NAAPS  may  contribute  to  the  uncertainty  of  the 
NAAPS  AOD  product.  However,  even  with  the  increase  in  uncertainty  at  high  AODs, 
Figure  39  reveals  that  NAAPS  provides  accurate  AOD  values,  on  an  annual  basis,  for  the 
Persian  Gulf  basin. 
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Figure  39.  Annual  scatter  plot  of  daily  AERONET  and  NAAPS  AOD  values  for  the 
Persian  Gulf  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one- 
to-one  line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size 
(N)  are  listed  in  the  upper  left  comer. 

Figure  40  contains  scatter  plots  of  AERONET  and  NAAPS  AOD  values  by 
season  for  the  2006-2007  timeframe.  The  correlation  is  consistent  across  seasons, 
indicating  that  the  accuracy  of  NAAPS  AOD  forecasts  varies  little  by  season. 
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Figure  40.  Persian  Gulf  Basin  AERONET  AOD  (0.50pm)  versus  NAAPS  AOD 

(0.55 jum)  scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c) 
Summer,  (d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line 
(blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are 

listed  in  the  upper  left  corner. 

In  conclusion,  the  results  for  the  Persian  Gulf  basin  NAAPS/AERONET 
comparison  are  very  promising.  Figures  38,  39,  and  40  suggest  that  the  NAAPS  with 
assimilated  MODIS  AOD  data  does  provide  representative  AOD  values  for  the  Persian 
Gulf  basin. 

2.  East  Asia 


a.  Korea 

The  NAAPS  AOD  grid  point  closest  to  Anmyon,  South  Korea  is  used  for 
the  Korea  basin  comparison.  Figure  41  is  the  two-year  NAAPS  0.55pm  AOD 
distribution  for  the  Korea  basin  for  2006-2007.  Half  of  the  AOD  values  lie  between  0.20 
and  0.45,  with  the  median  of  the  distribution  at  0.30.  In  75%  of  the  cases,  the  aerosol 
loading  may  be  light  enough  to  allow  for  scene  correction,  using  an  AOD  of  0.50  as  a 
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baseline.  In  comparison,  Korea  basin  AERONET  AOD  data  (Figure  16)  shows  that  (a) 
the  median  of  the  distribution  is  0.28  and  (b)  75%  of  the  cases  show  light  enough  aerosol 
loading  to  allow  for  scene  correction.  The  comparison  between  Figures  16  and  41 
suggests  that  for  the  two-year  period,  the  NAAPS  aerosol  loading  over  the  Korea  basin 
closely  matches  the  AERONET  aerosol  loading. 


Daily  Mean  NAAPS  AOD  (0.55  ^m) 


Figure  41 .  Two-year  distribution  of  0.55pm  NAAPS  AOD  values  for  Anmyon  in  the 
Korea  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 

Figure  42  contains  a  scatter  plot  of  AERONET  and  NAAPS  AOD  values 
for  the  Korea  basin  for  the  two-year  period  from  2006  through  2007.  The  scatter  is  not 
symmetric  around  the  one-to-one  line  for  AOD  values  above  0.80.  Figure  42  appears  to 
show  a  tendency  for  NAAPS  to  underestimate  high  AOD  events,  especially  when  the 
AOD  is  greater  than  0.80. 
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Figure  42.  Annual  scatter  plot  of  daily  AERONET  and  NAAPS  AOD  values  for  the 
Korea  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one-to-one 
line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N) 
are  listed  in  the  upper  left  corner. 

Figure  43  contains  scatter  plots  of  AERONET  and  NAAPS  AOD  values 
by  season  for  the  2006-2007  timeframe.  Correlation  coefficients  range  from  0.59  to  0.78. 
Figure  43  suggests  a  lower  correlation  between  NAAPS  and  AERONET  AOD  values  for 
the  Korea  basin  than  the  Persian  Gulf.  Additionally,  Fall  and  Winter  show  much  less 
scatter  than  Spring  and  Summer,  indicating  less  uncertainty  in  NAAPS  AOD  estimates  in 
Fall  and  Winter.  Finally,  the  correlation  is  highest  in  the  Winter. 

In  summary,  Figures  41,  42,  and  43  indicate  mixed  results  for  the  Korea 
basin.  The  NAAPS  distribution  over  the  two-year  period  (Figure  41)  compared  well  with 
the  AERONET  distribution  for  the  basin  (Figure  16).  The  point-to-point  daily  mean 
NAAPS  and  AERONET  comparison  indicates  a  tendency  for  NAAPS  to  underestimate 
AOD  on  days  with  AOD  values  greater  than  0.80. 
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Figure  43.  Korea  Basin  AERONET  AOD  (0.50pm)  versus  NAAPS  AOD  (0.55pm) 
scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c)  Summer, 
(d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line  (blue)  are 
denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are  listed  in  the 

upper  left  corner. 


b.  Japan 

The  NAAPS  AOD  grid  point  closest  to  each  of  the  two  Japan  stations, 
Shirahama  and  Osaka  is  used  in  the  Japan  basin  comparison.  Figure  44  is  the  two-year 
NAAPS  0.55pm  AOD  distribution  for  the  Japan  basin  for  2006-2007.  Half  of  the  AOD 
values  lie  between  0.15  and  0.30,  with  the  median  of  the  distribution  at  0.21.  In  95%  of 
the  cases,  the  aerosol  loading  may  be  light  enough  to  allow  for  scene  correction,  using  an 
AOD  value  of  0.50  as  a  baseline.  In  comparison,  Japan  basin  AERONET  AOD  data 
(Figure  17)  shows  that  (a)  the  median  of  the  distribution  is  also  0.21  and  (b)  87%  of  the 
cases  show  light  enough  aerosol  loading  to  allow  for  scene  correction.  The  comparison 
between  Figures  17  and  44  suggests  that  for  the  two-year  period,  the  NAAPS  estimates  of 
aerosol  loading  over  the  Japan  basin  closely  matches  the  AERONET  aerosol  loading. 
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Figure  44.  Two-year  distribution  of  0.55pm  NAAPS  AOD  values  for  Osaka  and 

Shirahama  in  the  Japan  basin,  for  2006-2007.  The  histogram  of  AOD  values 
(blue  bars)  and  cumulative  frequency  distribution  (red  curve)  corresponding  to 
the  histogram  are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer. 

The  vertical  lines  (black),  from  left  to  right,  denote  the  quartiles  of  the 
distribution,  25%  (dashed),  50%  (solid),  75%  (dashed). 

Figure  45  contains  a  scatter  plot  of  AERONET  and  NAAPS  AOD  values 
for  the  Japan  basin  for  the  two-year  period  from  2006  through  2007.  The  correlation  is 
0.73.  NAAPS  AOD  values  display  a  low  bias  at  higher  AOD  values  greater  than  0.40. 
Consequently,  NAAPS  AOD  values  for  the  Japan  basin  are  likely  to  be  low  when  aerosol 
loading  is  heavy.  Otherwise,  there  is  good  agreement  between  the  NAAPS  and 
AERONET  AOD  values. 
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Figure  45.  Annual  scatter  plot  of  daily  AERONET  and  NAAPS  AOD  values  for  the 
Japan  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one-to-one 
line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N) 
are  listed  in  the  upper  left  corner. 

Figure  46  contains  scatter  plots  of  AERONET  and  NAAPS  AOD  values 
by  season  for  2006-2007.  The  tendency  to  underestimate  higher  AOD  values,  seen  in 
Figure  45,  is  also  seen  in  the  seasonal  scatter  plots.  The  underestimation  is  most 
prominent  in  the  Spring  and  Summer  seasons  when  aerosol  loading  is  highest.  In  the  Fall 
season,  which  has  a  lighter  aerosol  load  than  Spring  and  Summer,  the  underestimation 
occurs  for  AOD  values  greater  than  0.25.  NAAPS  estimates  have  the  highest  correlation 
for  the  Winter  season.  Overall,  Figure  46  suggests  that  NAAPS  performs  best  over  the 
Japan  basin  when  aerosol  loading  is  low. 
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Figure  46.  Japan  Basin  AERONET  AOD  (0.50pm)  versus  NAAPS  AOD  (0.55pm) 
scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c)  Summer, 
(d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line  (blue)  are 
denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are  listed  in  the 

upper  left  corner. 

In  summary,  the  results  for  the  Japan  NAAPS/AERONET  comparison  are 
promising.  Figures  44,  45,  and  46  show  that  the  NAAPS  does  provide  representative 
AOD  values  for  the  Japan  basin  when  the  aerosol  loading  is  low. 

3.  West  Africa 

The  NAAPS  AOD  grid  point  closest  to  Ilorin,  Nigeria  is  used  for  the  West  Africa 
basin  comparison,  but  this  grid  point  is  over  water  approximately  350  km  due  south  of 
Ilorin.  Figure  47  is  the  two-year  NAAPS  0.55pm  AOD  distribution  for  the  West  Africa 
basin  for  2006-2007.  Half  of  the  AOD  values  lie  between  0.20  and  0.40,  with  the  median 
of  the  distribution  at  0.30.  In  82%  of  the  cases,  the  aerosol  loading  may  be  light  enough 
to  allow  for  scene  correction,  using  an  AOD  of  0.50  as  a  baseline.  In  comparison,  West 
Africa  basin  AERONET  AOD  data  (Figure  23)  shows  that  (a)  the  median  of  the 
distribution  is  0.70  and  (b)  only  35%  of  the  cases  show  light  enough  aerosol  loading  to 
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allow  for  scene  correction.  Therefore,  the  comparison  between  Figures  23  and  47 
suggest  that  on  an  annual  basis,  the  NAAPS  aerosol  loading  over  the  West  Africa  basin  is 
significantly  less  than  the  AERONET  aerosol  loading. 
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Figure  47.  Two-year  distribution  of  0.55|um  NAAPS  AOD  values  for  Ilorin  in  the  West 
Africa  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 

Figure  48  contains  a  scatter  plot  of  AERONET  and  NAAPS  AOD  values  for  the 
West  Africa  basin  for  the  two-year  period  from  2006  through  2007.  The  scatter  implies  a 
low  bias  in  the  NAAPS  AOD  values.  The  MODIS/AERONET  comparison  similarly 
reveals  a  low  bias  for  West  Africa.  Consequently,  the  NAAPS  AOD  estimates  for  the 
West  Africa  basin  are  likely  to  be  lower  than  AERONET  AOD  measurements. 
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Figure  48.  Annual  scatter  plot  of  daily  AERONET  and  NAAPS  AOD  values  for  the 
West  Africa  basin  for  2006-2007.  The  daily  AOD  values  (red  stars)  and  one- 
to-one  line  (blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size 
(N)  are  listed  in  the  upper  left  corner. 

Figure  49  contains  scatter  plots  of  AERONET  and  NAAPS  AOD  values  by 
season  for  2006-2007.  All  seasons  display  the  tendency  for  NAAPS  to  underestimate 
AOD  values.  The  correlation  is  not  consistent  across  seasons,  and  all  seasons  display  a 
low  bias. 

In  summary,  NAAPS  AOD  estimates  consistently  underestimate  aerosol  loading 
in  the  West  Africa  basin.  The  comparison  of  MODIS  and  AERONET  in  the  West  Africa 
basin  provided  similar  results.  With  both  the  NAAPS  and  MODIS  datasets,  the  over¬ 
water  grid  point  is  more  than  300  km  away  from  the  AERONET  station  and  dust  and 
pollution  sources.  It  is  reasonable  that  on  an  annual  basis  NAAPS  will  underestimate 
aerosol  loading. 
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Figure  49.  West  Africa  Basin  AERONET  AOD  (0.50pm)  versus  NAAPS  AOD 

(0.55  pm)  scatterplots  for  the  2006-2007  seasons:  (a)  Winter,  (b)  Spring,  (c) 
Summer,  (d)  Fall.  The  daily  AOD  values  (red  stars)  and  one-to-one  line 
(blue)  are  denoted.  The  correlation  coefficient  (R)  and  sample  size  (N)  are 
located  in  the  upper  left  comer. 


4.  Summary 

The  results  from  the  NAAPS/AERONET  comparison  indicate  that  the  NAAPS 
AOD  estimates  are  a  reliable  proxy  for  AERONET  AOD  data  over  the  Persian  Gulf  and 
Japan  basins,  especially  when  aerosol  loading  is  light.  The  comparison  results  for  the 
Korea  and  West  Africa  basins  are  less  conclusive,  due  to  low  sample  size  (Korea)  and 
distance  from  the  coast  (West  Africa).  Errors  in  the  modeled  wind  field,  aerosol  types, 
and  aerosol  sources  no  doubt  account  for  some  of  the  uncertainty  associated  with  the 
NAAPS  estimates.  Additionally,  there  is  uncertainty  in  the  MODIS  AOD  retrievals 
which  are  assimilated  into  NAAPS.  Limitations  of  the  MODIS  AOD  retrievals  are 
beyond  the  scope  of  this  thesis  and  are  not  investigated  here. 
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D.  RADIANCE-BASED  TARGET  DETECTION 


This  thesis  creates  a  sample  radiance-based  target  detection  matrix,  Table  5,  as 
described  in  Chapter  III,  Section  D.  The  matrix  simulates  the  fonn  of  information 
decision-makers  require  to  estimate  the  probability  that  a  satellite  sensor  can  distinguish  a 
ground-based  target  during  high  AOD  events.  The  matrix  is  constructed  using 
MODTRAN™  calculations,  and  displays  simulated  TO  A  radiance  values  by  AOD  and 
surface  reflectance  for  a  specific  wavelength  (0.63pm  in  this  study). 

To  validate  the  MODTRAN™  radiance  values  in  Table  5,  this  study  uses  the 
TOA  radiance  measurements  from  a  MODIS  image  at  three  locations  in  the  image.  For 
each  location,  the  AOD  and  surface  reflectance  are  estimated.  Figure  50  is  a  MODIS 
image  for  0630  UTC  14  May  2004,  which  contains  a  large  dust  plume  extending  from 
southern  Iraq  down  to  western  portions  of  United  Arab  Emirates.  Point  1  in  Figure  50 
represents  a  location  with  a  low  amount  of  dust.  Point  2  is  a  location  with  a  moderate 
amount  of  dust,  and  Point  3  is  a  location  with  a  high  amount  of  dust.  Using  the 
TeraVision  satellite  image  display  software,  TOA  radiance  values  are  measured  at  each 
point  for  the  following  wavelengths:  0.44pm,  0.55pm,  0.65pm,  0.86pm,  and  1.60pm. 
The  measured  MODIS  TOA  radiation  values  at  Points  1,  2,  and  3  are  listed  in  Table  6. 


66 


Table  5.  MODTRAN™-Simulated  TOA  Radiance  (Wm'jum'sr  ' )  at  0.63|um  as  a 
function  of  AOD  (red)  and  Surface  Reflectance  (%,  blue). 
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Figure  50.  14  May  2004  0630  UTC  MODIS  image  for  the  Persian  Gulf.  Point  1  is  a 

location  with  a  low  amount  of  dust,  Point  2:  moderate  dust,  and  Point  3:  high 

dust. 
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Table  6.  MODIS  Top-of-the -Atmosphere  Radiance  Values  for  Points  1,  2,  and  3  in 

Figure  50.  Radiance  units  are  Wm  'jum  'sr  ' . 


Wavelength 

(pm) 

Point  1 
Radiance 
(Low  dust) 

Point  2 

Radiance 
(Moderate  dust) 

Point  3 

Radiance 
(High  dust) 

0.44 

80.40 

98.88 

140.35 

0.55 

43.91 

69.01 

141.92 

0.65 

27.14 

49.30 

126.91 

0.86 

12.85 

25.81 

77.99 

1.60 

2.39 

4.45 

15.30 

Comparison  of  the  radiance  values  for  Points  1,  2,  and  3  reveals  that  as 
atmospheric  aerosol  (dust)  content  increases,  so  does  the  TOA  radiance  received  by  the 
satellite.  The  dust  in  the  atmosphere  scatters  incoming  solar  radiation,  much  of  which  is 
reflected  back  to  the  satellite.  Therefore,  the  TOA  radiance  values  increase  as 
atmospheric  aerosol  content  increases. 

Figure  51  provides  AOD  estimates  at  0.63pm  for  the  Persian  Gulf  using  NOAA- 
17  imagery  at  0730  UTC  14  May  04  (Durkee  et  al.  1991;  NRL-Monterey  2009).  The 
AOD  values  in  Figure  51  are  an  estimate  of  the  AOD  values  in  the  MODIS  image  in 
Figure  50.  The  image  contains  AOD  values  for  the  eastern  half  of  the  Persian  Gulf.  No 
AOD  values  for  the  western  Persian  Gulf  are  available  because  the  dust  is  so  thick  that 
the  AOD  retrieval  algorithm  interprets  the  dust  as  cloud  cover  (Martins  et  al.  2002). 
The  color-coded  legend  on  the  left  side  of  the  image  lists  the  AOD  values.  The  AOD 
values  at  Point  1,  2,  and  3  are  respectively:  0.5,  1.5,  and  2.4. 
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Figure  51.  14  May  2004  0730  UTC  NOAA-17  image,  0.63|um,  denoting  retrieved 

AODs  for  the  Persian  Gulf  basin  (NRL-Monterey,  2009).  Points  1,  2,  and  3 
are  the  same  locations  denoted  in  Figure  50. 


A  representative  ocean  surface  reflectance  for  Points  1-3  is  five  percent  (ASTER 
Spectral  Library  2009).  Using  the  five  percent  surface  reflectance  column  in  Table  5  and 
the  AOD  values  corresponding  to  Points  1-3,  the  MODTRAN™-estimated  TO  A  radiance 
values  are  extracted  from  Table  5.  The  MODIS  TOA  radiance  measurements  at  0.65  pm 
and  the  MODTRAN™  TOA  radiance  estimates  at  0.63  pm  are  displayed  in  Table  7. 


Table  7.  TOA  Radiance  ( Wm  'jum  'sr  '  )  from  the  MODIS  image  (Figure  50)  and 
MODTRAN™  simulation  (Table  5)  and  the  NOAA-17  AOD  estimates  from 
Figure  5 1 . 


Point 

NOAA-17  AOD 
(0.63  pm) 

MODIS 
(0.65  pm) 

MODTRAN™ 
(0.63  pm) 

1 

0.50 

27.14 

27.10 

2 

1.50 

49.30 

38.50 

3 

2.40 

126.91 

47.50 
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At  Point  1,  there  is  little  difference  between  the  measured  and  simulated  TO  A 
radiance  values.  At  Point  2,  with  moderate  dust,  the  MODIS  TOA  radiance  is  higher 
than  the  simulated  radiance.  At  Point  3,  with  high  dust  loading,  significantly  more 
radiance  is  detected  by  the  MODIS  sensor  than  MODTRAN™  simulates.  As  AOD 
increases  above  approximately  0.5,  MODTRAN™  may  not  fully  account  for  the  amount 
of  multiple  scatter  interactions  with  the  aerosols,  which  would  reduce  the  simulated  TOA 
radiance.  Consequently,  MODTRAN™  estimates  of  TOA  radiance  are  currently 
insufficient  to  produce  a  matrix  for  operational  use. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FUTURE 

RESEARCH 


A.  CONCLUSIONS 

The  MODIS  data  used  in  this  study  is  reliable  in  retrieving  AOD  values  at 
0.55pm  for  the  Persian  Gulf  and  Japan  basins,  when  compared  to  AERONET  data  at 
0.50jiun.  Consequently,  data  from  this  MODIS  retrieval  may  be  used  with  confidence 
over  the  Persian  Gulf  and  Japan.  For  the  Korea  basin,  the  MODIS  AOD  retrieval  at 
0.55pm  may  provide  close  to  ground  truth  AOD  values.  Larger  AERONET  and  MODIS 
AOD  sample  sizes  will  help  reduce  uncertainty  and  increase  the  reliability  of  MODIS 
AOD  measurements  over  the  Korea  basin.  As  for  West  Africa,  MODIS-retrieved  AOD 
values  consistently  underestimate  aerosol  loading.  As  such,  over-water  MODIS-retrieved 
AOD  values  are  not  a  reliable  proxy  for  AERONET  AOD  values  over  continental  West 
Africa. 

Likewise,  the  NAAPS  with  NAVDAS-AOD  model  used  in  this  study  is  reliable 
in  forecasting  AOD  values  at  0.55pm  for  the  Persian  Gulf  and  Japan  basins,  when 
compared  to  AERONET  data  at  0.50pm.  Consequently,  data  from  this  model  may  be 
used  with  confidence  over  the  Persian  Gulf  and  Japan.  The  accuracy  of  NAAPS 
forecasts  for  Korea  and  West  Africa  is  less  certain.  NAAPS  estimates  show  some  skill 
for  Korea,  but  larger  sample  sizes  should  improve  accuracy  over  Korea.  The  low  bias 
over  West  Africa  is  most  likely  due  to  the  distance  between  the  NAAPS  AOD  grid  point 
and  the  AERONET  station. 

This  thesis  is  unable  to  validate  the  MODTRAN™  TOA  radiance  calculation 
because  of  the  large  disparity  between  the  MODTRAN™  radiances  and  MODIS 
radiances.  Consequently,  MODTRAN™  estimates  of  TOA  radiance  are  currently 
insufficient  to  produce  a  matrix  for  operational  use. 
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B.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

This  study  uses  NAAPS  and  MODIS  data  at  only  one  wavelength.  More  research 
needs  to  be  done  to  detennine  how  well  each  system  characterizes  aerosols  at  other 
wavelengths  and  thereby  detennine  the  impacts  of  different  particle  sizes. 

Additionally,  lack  of  data  was  a  problem  over  Korea.  In  the  future,  data  from 
eastern  China  needs  to  be  used  to  help  account  for  conditions  on  the  east  coast  of  the 
Yellow  Sea.  Also,  if  new  AERONET  stations  happen  to  be  set  up  along  the  Atlantic 
Coast  in  West  Africa,  their  data  needs  to  be  included  in  future  studies  for  this  region. 
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APPENDIX  A:  AOD  VARIATION  BY  WAVELENGTH 


Figure  52.  Two-year  distribution  of  AERONET  AOD  (0.44|um)  values  for  Solar  Village 
for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and  cumulative 
frequency  distribution  (red  curve)  corresponding  to  the  histogram  are  denoted. 
The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical  lines  (black), 
from  left  to  right,  denote  the  quartiles  of  the  distribution,  25%  (dashed),  50% 

(solid),  75%  (dashed). 
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100% 


AERONET  AOD  (0.50  pm) 

Figure  53.  Same  as  Figure  52,  but  for  0.50|um. 


AERONET  AOD  (0.675pm) 


Figure  54.  Same  as  Figure  52,  but  for  0.675|um. 
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Cumulative  Frequency  (percent) 


Cumulative  Probability  (%) 


AERONET  AOD  (0.870|im) 


Figure  55.  Same  as  Figure  52,  but  for  0.870|um. 


Figure  56.  Solar  Village  AERONET  AOD  Cumulative  Frequency  Curves  for  0.44|um, 

0.50jum,  0.675)iim  and  0.870|um  for  2006-2007. 
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APPENDIX  B:  AERONET  CLIMATOLOGY 


A.  PERSIAN  GULF 


AERONET  AOD  (0.50  (im) 
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o 


Figure  57.  One -year  distribution  of  0.50  |um  AERONET  AOD  values  for  the  Persian 
Gulf  basin  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  58.  Same  as  Figure  57,  but  for  2007. 


AERONET  AOD  (0-50  p/n) 

Figure  59.  Two-year  Winter  cumulative  frequency  curves  for  the  Persian  Gulf  Basin  for 
2006-2007.  The  cumulative  frequency  curve  for  each  individual  station  is 
compared  with  the  cumulative  frequency  curve  computed  from  the  histogram 
distribution  for  the  combined  dataset  of  the  three  stations  (red  curve). 
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Figure  60.  Same  as  Figure  59,  but  for  Spring. 


Figure  6 1 .  Same  as  Figure  59,  but  for  Summer. 
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Figure  62.  Same  as  Figure  59,  but  for  Fall. 
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Two-year  Winter  cumulative  frequency  curves  for  Korea  (blue)  and  Japan 

(red)  for  2006-2007. 
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Figure  64.  Same  as  Figure  63,  but  for  Spring. 
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Figure  65.  Same  as  Figure  63,  but  for  Summer. 


Figure  66.  Same  as  Figure  63,  but  for  Fall. 
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Cumulative  Probability  (%)  Cumulative  Probability (%) 


Two-year  Winter  cumulative  frequency  curves  for  the  Japan  Basin  for  2006- 
2007.  The  cumulative  frequency  curve  for  each  individual  station  is 
compared  with  the  cumulative  frequency  curve  computed  from  the  histogram 
distribution  for  the  combined  dataset  of  the  two  stations  (red  curve).  Here,  no 
data  is  available  from  Osaka,  so  only  Shirahama  is  shown. 
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Figure  68.  Same  as  Figure  67,  but  for  Spring. 
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Figure  69.  Same  as  Figure  67,  but  for  Summer. 


Figure  70.  Same  as  Figure  67,  but  for  Fall. 
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One-year  distribution  of  0.50  (am  AERONET  AOD  values  for  the  Korea 
basin  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and  cumulative 
frequency  distribution  (red  curve)  corresponding  to  the  histogram  are  denoted. 
The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical  lines  (black), 
from  left  to  right,  denote  the  quartiles  of  the  distribution,  25%  (dashed),  50% 

(solid),  75%  (dashed). 


AERONET  AOD  (0.50  (im) 


Figure  72.  Same  as  Figure  71,  but  for  2007. 
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AERONET  AOD  (0.50  |im) 


One-year  distribution  of  0.50  pm  AERONET  AOD  values  for  the  Japan 
basin  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and  cumulative 
frequency  distribution  (red  curve)  corresponding  to  the  histogram  are  denoted. 
The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical  lines  (black), 
from  left  to  right,  denote  the  quartiles  of  the  distribution,  25%  (dashed),  50% 

(solid),  75%  (dashed). 
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Figure  74.  Same  as  Figure  73,  but  for  2007. 
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Figure  75.  One-year  distribution  of  0.44|um  and  0.50|um  AERONET  AOD  values  for  the 
West  Africa  basin  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  77.  Two-year  Winter  cumulative  frequency  curves  for  the  West  Africa  Basin  for 
2006-2007.  The  cumulative  frequency  curve  for  each  individual  station  is 
compared  with  the  cumulative  frequency  curve  computed  from  the  histogram 
distribution  for  the  combined  dataset  of  the  two  stations  (red  curve). 
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Figure  78.  Same  as  Figure  77,  but  for  Spring. 


Figure  79.  Same  as  Figure  77,  but  for  Summer. 
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Figure  80.  Same  as  Figure  77,  but  for  Fall. 


AERONET  AOD  (0.44  ^m) 


Figure  8 1 .  Winter  distribution  of  AERONET  AOD  (0.44|um)  values  for  station  Djougou 

in  the  West  Africa  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue 
bars)  and  cumulative  frequency  distribution  (red  curve)  corresponding  to  the 
histogram  are  denoted,  with  AOD  values  ranging  from  0  to  4.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  82.  Same  as  Figure  8 1 ,  but  for  Spring. 


Figure  83.  Same  as  Figure  81,  but  for  Summer. 
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Figure  84.  Same  as  Figure  81,  but  for  Fall. 
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APPENDIX  C:  MODIS/AERONET  COMPARISON 


A.  PERSIAN  GULF 
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Figure  85.  Two-year  distribution  of  0.55  |um  MODIS  AOD  values  for  stations  Solar 
Dhabi  and  Dhadnah  in  the  Persian  Gulf  basin,  for  2006-2007.  The  histogram 
of  AOD  values  (blue  bars)  and  cumulative  frequency  distribution  (red  curve) 
corresponding  to  the  histogram  are  denoted.  The  sample  size  is  listed  in  the 
upper  right  corner.  The  vertical  lines  (black),  from  left  to  right,  denote  the 
quartiles  of  the  distribution,  25%  (dashed),  50%  (solid),  75%  (dashed). 
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Figure  86.  Seasonal  cumulative  frequency  curves  for  the  Persian  Gulf  basin  MODIS 
(red)  and  AERONET  (blue)  AOD  distributions  for  the  2006-2007  seasons: 

(a)  Winter,  (b)  Spring,  (c)  Summer,  (d)  Fall.  The  y-axis  indicates  the 
cumulative  probability  of  occurrence  of  an  AOD.  The  x-axis  indicates  the 

AOD  values. 


94 


AOD  -  Count  per  Bin 


B. 


EAST  ASIA 


1.  Korea 


MODIS  AOD  (0.55  (im) 


E 

o 


Figure  87.  Two-year  distribution  of  0.55  jam  MODIS  AOD  values  for  Anmyon  in  the 
Korea  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  88.  Seasonal  cumulative  frequency  curves  for  the  Korea  basin  MODIS  (red)  and 
AERONET  (blue)  AOD  distributions  for  the  2006-2007  seasons:  (a)  Winter, 
(b)  Spring,  (c)  Summer,  (d)  Fall.  The  y-axis  indicates  the  cumulative 
probability  of  occurrence  of  an  AOD.  The  x-axis  indicates  the  AOD  values. 
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Figure  89.  Two-year  distribution  of  0.55  pm  MODIS  AOD  values  for  Shirahama  in  the 
Japan  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  90.  Seasonal  cumulative  frequency  curves  for  the  Japan  basin  MODIS  (red)  and 
AERONET  (blue)  AOD  distributions  for  the  2006-2007  seasons:  (a)  Winter, 
(b)  Spring,  (c)  Summer,  (d)  Fall.  The  y-axis  indicates  the  cumulative 
probability  of  occurrence  of  an  AOD.  The  x-axis  indicates  the  AOD  values. 
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Figure  91.  Two-year  distribution  of  0.55  jam  MODIS  AOD  values  for  Ilorin  in  the  West 
Africa  basin,  for  2006-2007.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  92.  Seasonal  cumulative  frequency  curves  for  the  West  Africa  basin  MODIS 
(red)  and  AERONET  (blue)  AOD  distributions  for  the  2006-2007  seasons: 

(a)  Winter,  (b)  Spring,  (c)  Summer,  (d)  Fall.  The  y-axis  indicates  the 
cumulative  probability  of  occurrence  of  an  AOD.  The  x-axis  indicates  the 

AOD  values. 
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APPENDIX  D:  NAAPS/AERONET  COMPARISON 


A.  PERSIAN  GULF 


Daily  Mean  NAAPS  AOD  (0.55  jim) 


Figure  93.  One-year  distribution  of  0.55jum  NAAPS  AOD  values  for  stations  Solar 
Village,  Dhabi  and  Dhadnah  in  the  Persian  Gulf  basin,  for  2006.  The 
histogram  of  AOD  values  (blue  bars)  and  cumulative  frequency  distribution 
(red  curve)  corresponding  to  the  histogram  are  denoted.  The  sample  size  is 
listed  in  the  upper  right  comer.  The  vertical  lines  (black),  from  left  to  right, 
denote  the  quartiles  of  the  distribution,  25%  (dashed),  50%  (solid),  75% 

(dashed). 
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Figure  94.  Same  as  Figure  93,  but  for  2007. 


Daily  Mean  NAAPS  AOD  (0.55  ^m) 

Figure  95.  Seasonal  distribution  of  0.55|um  NAAPS  AOD  values  for  the  Persian  Gulf 
basin,  for  2006-2007.  The  seasons  and  their  respective  median  values  of  the 
distribution  are  located  in  the  upper  right  corner.  The  y-axis  indicates  the 
number  of  AOD  samples  per  bin  from  0  to  30,  with  the  x-axis  denoting  AOD 

values  at  0.55pm. 
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Figure  96.  One -year  distribution  of  0.5  5  (am  NAAPS  AOD  values  for  Anmyon  in  the 

Korea  basin,  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  97.  Same  as  Figure  96,  but  for  2007. 
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Figure  98.  Seasonal  distribution  of  0.55|um  NAAPS  AOD  values  for  the  Korea  basin, 
for  2006-2007.  The  seasons  and  their  respective  median  values  of  the 
distribution  are  located  in  the  upper  right  corner.  The  y-axis  indicates  the 
number  of  AOD  samples  per  bin  from  0  to  30,  with  the  x-axis  denoting  AOD 

values  at  0.5 5 pm. 
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Daily  Mean  NAAPS  AOD  (0.55  jim) 


Figure  99.  One-year  distribution  of  0.55pm  NAAPS  AOD  values  for  Osaka  and 

Shirahama  in  the  Japan  basin,  for  2006.  The  histogram  of  AOD  values  (blue 
bars)  and  cumulative  frequency  distribution  (red  curve)  corresponding  to  the 
histogram  are  denoted.  The  sample  size  is  listed  in  the  upper  right  corner.  The 
vertical  lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution, 
25%  (dashed),  50%  (solid),  75%  (dashed). 
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Figure  100.  Same  as  Figure  99,  but  for  2007. 


Daily  Mean  NAAPS  AOD  (0.55  jim) 

Figure  101.  Seasonal  distribution  of  0.55pm  NAAPS  AOD  values  for  the  Japan  basin,  for 
2006-2007.  The  seasons  and  their  respective  median  AOD  values  of  the 
distribution  are  located  in  the  upper  right  corner.  The  y-axis  indicates  the 
number  of  AOD  samples  per  bin  from  0  to  30,  with  the  x-axis  denoting  AOD 

values  at  0.55pm. 
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Daily  Mean  NAAPS  AOD  (0.55  yim) 


Figure  102.  One-year  distribution  of  0.55|um  NAAPS  AOD  values  for  Ilorin  in  the  West 
Africa  basin,  for  2006.  The  histogram  of  AOD  values  (blue  bars)  and 
cumulative  frequency  distribution  (red  curve)  corresponding  to  the  histogram 
are  denoted.  The  sample  size  is  listed  in  the  upper  right  comer.  The  vertical 
lines  (black),  from  left  to  right,  denote  the  quartiles  of  the  distribution,  25% 
(dashed),  50%  (solid),  75%  (dashed). 
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Figure  103.  Same  as  Figure  102,  but  for  2007. 


Daily  Mean  NAAPS  AOD  (0.55  jim) 

Figure  104.  Seasonal  distribution  of  0.5 5 (am  NAAPS  AOD  values  for  the  West  Africa, 
for  2006-2007.  The  seasons  and  their  respective  median  values  of  the 
distribution  are  located  in  the  upper  right  corner.  The  y-axis  indicates  the 
number  of  AOD  samples  per  bin  from  0  to  30,  with  the  x-axis  denoting  AOD 

values  at  0.55pm. 
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